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INTRODUCTION 


The  interrogation  of  tissue  obtained  from  patients  with  castration-resistant  prostate  cancer  (CRPC)  reveal  molecular 
changes  associated  with  this  disease  state;  however,  our  understanding  of  signaling  pathways  functionally  leading  to 
castration  resistance  remains  incomplete.  Here  we  have  performed  a  high  throughput,  in  vivo  genetic  screen  to  identify 
kinases  that  permit  androgen-dependent  transformed  prostate  epithelial  cells  (LHSR-AR  cells)  to  form  tumors  in  female 
animals.  In  addition  to  known  prostate  cancer  oncogenes  and  mediators  of  androgen  independence  (mutated  KRAS, 
constitutively  active  MEK,  RAF  I,  ERBB2,  AKTI,  PIMI  and  PIM2),  overexpression  of  the  Never  In  Mitosis  A  (NIMA) 
related  kinase  6  (NEK6)  reproducibly  yielded  androgen-independent  tumors.  NEK6  is  overexpressed  in  prostate  cancer 
cell  lines  compared  to  their  normal  counterparts  and  is  overexpressed  in  a  subset  of  human  prostate  cancers.  Expression 
of  NEK6  confers  castration  resistance  to  established  tumors  in  male  mice,  and  suppressing  NEK6  expression  restores 
sensitivity  to  castration.  Castration-resistant  tumors  generated  through  NEK6  overexpression  are  predominantly  squamous 
in  histology  and  do  not  express  androgen  receptor  (AR),  and  NEK6  does  not  activate  AR  signaling.  Phosphoproteome 
and  interactome  analysis  reveals  the  transcription  factors  FOXJ2  and  NCOA5,  as  well  as  the  kinases  CKla  and  YESl  to 
be  novel  substrates.  The  gene  expression  profile  mediated  by  NEK6  overexpression  in  tumors  from  castrated  mice 
demonstrates  elements  of  both  differentiation  and  immune  signaling,  and  NEK6  maintains  the  expression  of  a  set  of 
genes,  including  genes  involved  in  interferon  signaling,  decreased  with  castration  in  this  model.  Phosphomimic  forms  of 
FOXJ2,  YESl,  and  CKla  (but  not  NCOA5)  recapitulate  elements  of  the  NEK6  gene  expression  signature,  and  the 
phosphomimic  CKla  signature  in  particular  correlates  highly  with  the  NEK6  signature  with  regards  to  expression  of 
markers  of  squamous  differentiation,  such  as  KRT13.  These  studies  reveal  NEK6/CKla  signaling  as  one  mechanism  of 
resistance  in  androgen  pathway  independent  prostate  cancer  (APIPC). 

BODY 

Background 

Prostate  cancer  is  the  second  most  common  cause  of  cancer  death  in  men,  and  the  majority  of  these  deaths  occur 
in  patients  with  metastatic  castration-resistant  prostate  cancer  (CRPC).  Clinical  responses  to  agents  that  decrease 
circulating  androgens  to  below  castrate  levels  (abiraterone  -  Attard  et  al,  2009)  and  new  potent  antagonists  of  the 
androgen  receptor  (enzalutamide  -  Scher  et  aL,  2012)  demonstrate  that  the  androgen  receptor  (AR)  signaling  pathway 
remains  critical  in  CRPC.  Recent  reports  from  studies  of  metastatic  tissue  in  patients  who  have  progressed  on  these 
therapies  reveal  that  tumors  can  become  resistant  either  through  persistent  activation  of  the  AR  pathway,  or  by 
progression  to  a  state  described  as  androgen  pathway  independent  prostate  cancer  (APIPC)  (Nelson,  2012),  with  or 
without  evidence  for  neuroendocrine  transdifferentiation. 

Constitutive  activation  of  kinases  such  as  ERBB2,  MAPK,  PBK/Akt,  and  Src  (Edwards  and  Bartlett,  2005b)  has 
been  implicated  in  mediating  resistance  to  hormonal  therapies  through  both  AR-dependent  (Lin  et  ai,  2001;  Yeh  et  al., 
1999;  Guo  et  ai,  2006),  and  AR-independent  mechanisms  (Pienta  and  Bradley,  2006).  A  number  of  recent  studies  have 
demonstrated  the  complementarity  and  crosstalk  between  kinase  signaling  pathways  and  AR  signaling.  Mendiratta,  et  al. 
(2009)  developed  a  gene  expression  signature  of  AR  signaling,  and  found  that  decreased  predicted  AR  activity  in  patient 
samples  correlated  with  increased  predicted  Src  activity.  Carver,  et  al.  (201 1)  and  Mulholland,  et  al.  (201 1)  demonstrated 
cross-regulation  through  reciprocal  feedback  between  the  AR  and  PTEN/PI3K/MTOR  signaling  pathways.  However, 
inhibitors  of  many  of  these  kinases  have  failed  to  demonstrate  significant  clinical  benefit  in  unselected  patient 
populations,  such  as  in  trials  of  the  ERBB2  inhibitor  lapatinib  (Whang  et  al.,  2011),  the  MTOR  inhibitor  everolimus 
(Nakabayashi  et  al.,  2012),  and  the  SRC  inhibitor  dasatinib  (Araujo  et  al.,  2013).  An  agent  targeting  MET  and  VEGFR2, 
cabozantinib,  has  demonstrated  interesting  clinical  activity  in  metastatic  CRPC  (Smith  et  al.,  2014),  but  has  not 
demonstrated  improvements  in  overall  survival  or  bone  pain  in  randomized  phase  III  trials. 

Tumor  samples  from  patients  with  CRPC  have  increased  levels  of  tyrosine  phosphorylation  as  compared  to 
treatment  naive  prostate  cancer  (Drake,  et  al.  2012;  Drake  et  al.,  2013),  and  analysis  of  serine/threonine  and  tyrosine 
phosphorylation  events  in  metastatic  CRPC  suggests  particular  signaling  networks  involved  in  this  transition  (Yazdi,  et  al. 
2014).  However,  there  is  limited  evidence  for  activating  kinase  point  mutations  in  CRPC,  suggesting  that  kinase  pathways 
are  activated  by  other  (structural  genetic,  epigenetic,  microenvironmental)  mechanisms.  Kinase  signaling  pathways 
whose  activation  could  lead  to  castration  resistance  in  patients  have  not  been  comprehensively  catalogued,  so  we  have 
performed  an  in  vivo  functional  genomic  screen  to  identify  novel  pathways  that  may  be  involved.  These  studies  would 
help  identify  novel  targets  for  therapeutic  intervention  (either  the  kinases  themselves,  or  upstream/downstream  mediators) 
and  help  with  biomarker  development  for  stratifying  patients  for  likelihood  of  response  to  therapy. 
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We  hypothesize  that  novel  signaling  pathways  may  be  involved  in  conferring  castrate  resistance,  and  given  that 
kinases  usually  act  as  transducers  of  growth  and  proliferation  signals,  we  hypothesize  more  specifically  that 
activated/amplified  kinases  play  a  role  in  the  development  of  castration  resistance.  We  have  thus  performed  an  in  vivo 
functional  genomic  screen  to  identify  novel  pathways  that  may  be  involved  and  likely  serve  as  therapeutic  targets  in  these 
patients. 

Previous  work  in  our  laboratory  (Berger  et  aL,  2004)  demonstrated  that  primary  prostate  epithelial  cells  (PrECs) 
that  are  rendered  tumorigenic  by  the  expression  of  the  SV40  large  T  and  small  t  antigens,  the  catalytic  subunit  of 
telomerase  (hTERT),  H-Ras,  and  the  androgen  receptor  (LHSR-AR  cells)  form  well-differentiated  tumors  in  mice.  These 
tumors  are  androgen-dependent  and  are  thus  unable  to  grow  in  female  or  castrated  mice.  I  have  performed  a  high 
throughput,  in  vivo  genetic  screen  to  identify  kinases  that  permitted  these  cells  to  form  tumors  in  female  animals.  A 
lentivirally  delivered  kinase  ORE  library  encompassing  601  kinases  and  other  oncogenes  was  introduced  into  these  cells 
in  pools  of  9-10,  and  I  identified  fourteen  ORE  integrants  that  allowed  for  the  androgen- independent  development  of 
subcutaneous  tumors  by  PCR  using  vector  specific  primers  (Table  1).  Using  the  same  ORE  library,  I  performed  an  in 
vitro  screen  for  genes  conferring  androgen-independent  proliferation  to  androgen  dependent  LNCaP  cells  and  identified 
13  genes  that  significantly  (>2  standard  deviations  from  median)  increased  proliferation  in  androgen-deprived  conditions 
(Table  1). 

The  24  total  candidates  identified  from  both  screens  were  introduced  into  LHSR-AR  cells  individually  and 
injected  into  6  subcutaneous  sites  of  female  B ALB/C  nude  mice  for  validation  of  androgen-independent  tumor  formation. 
Among  the  candidates  that  reproducibly  yield  androgen- independent  tumors  are  mutated  KRAS;  RAFl,  which  is 
recurrently  translocated  (Palanisamy  et  al.  2010)  and  amplified  (Taylor  et  ai,  2010)  in  advanced  prostate  cancer;  ERBB2, 
AKTl,  and  constitutively  active  MEKl,  which  have  been  implicated  in  androgen  independence  (Edwards  and  Bartlett, 
2005a);  and  PIMl  and  PIM2,  which  have  previously  been  demonstrated  to  be  important  oncogenes  in  prostate  cancer 
(Brault  et  al.,  2010).  Among  the  strongest  candidates  identified  to  confer  androgen  independence  in  this  assay  are  the 
Never  In  Mitosis  A  (NIMA)  related  kinase  6  (NEK6),  and  nemo-like  kinase  (NEK). 

Aim  1.  Elucidating  the  role  of  NIMA-r  elated  kinase  6  (NEK6)  and  nemo-like  kinase  (NLK)  as  mediators  of  castrate- 
resistant  prostate  cancer  and  assessing  their  potential  as  therapeutic  targets. 

Milestone  1:  Determine  whether  NEK6  and  NLK  can  confer  castrate  resistance  to  established  tumors,  whether  kinase 
activity  is  required,  and  whether  their  continued  expression  is  required. 

In  addition  to  conferring  tumor  formation  in  female  mice,  overexpression  of  NEK6  in  LHSR-AR  cells  also  lead  to 
tumor  formation  in  castrated  mice, 
which  lack  circulating  androgens 
since  mice  do  not  synthesize 
androgens  from  their  adrenal 
glands  (Figure  lA).  The  tumors  in 
castrated  mice  are  generally 
smaller  and  slower  growing  than 
those  in  female  mice,  suggesting 
that  signaling  mediated  by 
androgens  play  a  role  in  the 
growth  of  the  tumors  even  though 
they  are  not  essential  for  tumor 
formation.  In  addition,  we  tested 
whether  NEK6  overexpression 
can  lead  to  androgen-independent 
tumor  formation  in  a  different 
genetic  context,  that  is  in 
immortalized  PrECs  transformed 
by  overexpression  of  the  MYC 
oncogene  and  expression  of  an 
active  form  of  the  pi  10a  subunit 
of  PI3  kinase  (rather  than  with 


Table  1.  Results  of  in  vivo  and  in  vitro  screens  for  genes  conferring  androgen 
independence,  and  of  in  vivo  validation. 

In  vivo  LHSR-AR  screen 

In  vitro  LNCaP  screen 

In  vivo  validation 

ORF 

# 

tumors 

ORF 

Fold 

proliferation 
(x  median)^ 

ORF 

# 

tumors 

KRASV12+MEKDD‘ 

^  3/3 

NLK 

2.09 

KRASV12 

6/6 

ERBB2 

3/3 

CDK6 

1.90 

ERBB2 

5/6 

NEK6 

3/3 

PIMl 

1.72 

NEK6 

5/6 

RAFl 

3/3 

CDK4 

1.55 

NLK 

4/6 

AKTl 

1/3 

P1M2 

1.45 

MEKDD 

4/6 

BRD3+NEK8" 

1/3 

STK40 

1.44 

AKTl 

3/6 

P1M2 

1/3 

RPS6KA2 

1.40 

RAFl 

3/6 

GK 

1/3 

AGK 

1.40 

PIMl 

3/6 

PFKP 

1/3 

TGFBRl 

1.40 

P1M2 

1/6 

PRKG2 

1/3 

DAPK3 

1.39 

others 

0/6 

TGFBR2 

1/3 

LOC389599 

1.38 

PlMl 

1/3 

NEKS 

1.37 

AKTl 

1.37 

^  Two  ORF  integrants  were  amp] 

iified  from  these  tumors 

^  For  reference,  the  synthetic  androgen  R1881  leads  to  median  2.06  fold 
proliferation  ^ 


small  t  antigen  and  active  HRAS),  termed  LHMK-AR  cells.  Indeed,  NEK6  overexpression  led  to  tumor  formation  in  this 
other  genetic  context  (Figure  IB);  however  given  that  the  tumors  were  smaller  and  somewhat  less  robust  than  those 
formed  from  LHSR-AR  cells,  functional  studies  performed  below  focused  on  tumors  formed  from  the  latter  cells. 

.  As  reported  in  previous  Annual  Summary  Reports,  the  kinase  activity  of  NEK6  is  required  for  castration- 
resistant  tumor  formation,  as  a  kinase  dead  version  with  mutations  of  lysine  residues  at  amino  acids  74  and  75  to 
methionine  (K74M/K75M)  did  not  confer  this  phenotype.  To  assess  whether  NEK6  could  serve  as  a  therapeutic  target  in 
established  tumors  where  its  activity  is  increased,  we  implanted  cells  with  doxycycline-inducible  expression  of  NEK6 
(Figure  1C)  in  male  mice  with  a  testosterone  pellet,  and  tumors  were  allowed  to  form  in  the  presence  of  doxycycline  in 
the  diet.  After  35  days,  the  mice  were  castrated,  testosterone  pellets  was  removed,  and  doxycycline  was  either  continued 
or  withdrawn  (5  mice  =15  tumors  per  group;  mice  that  died  perioperatively  were  excluded  from  the  analysis).  We  found 
that  at  day  30  after  castration  NEK6  confers  resistance  to  castration  compared  to  the  parental  cells  when  its  expression  is 
maintained  with  doxycycline  (p=0.001)  but  sensitivity  to  castration  is  restored  when  doxycycline  is  withdrawn  (p=0.049) 
(Figure  1C). 

The  tumors  formed  due  to  NEK6  overexpression  have  regions  of  nuclear  AR  expression  but  the  majority  of  these 
tumors  lacked  AR  expression  (Figure  2A).  Hematoxylin  and  eosin  (H&E)  staining  demonstrates  squamous  differentiation 
in  these  tumors  with  the  more  mature  differentiated  regions  demonstrating  keratin  deposition  and  AR  loss.  Androgen- 
independent  tumors  derived  from  expression  of  the  other  kinases  identified  in  the  screen  also  have  regions  of  AR- 
positivity  and  negativity  that  vary  in  proportion  and  intensity  (Supplemental  Table  1).  NEK6  overexpressing  tumors 
established  in  male  mice  demonstrate  strong  nuclear  AR  staining  with  no  squamous  differentiation  (Figure  2B,  left 
panels).  However,  after  castration,  these  tumors  lose  nuclear  AR  over  time  and  develop  more  keratinization  between  nests 
of  tumor  cells  (Figure  2B,  bottom  panels).  This  observation  suggests  that  NEK6-mediated  castration  resistance  does  not 
require  AR  activity,  and  that  transdifferentiation  to  a  squamous  phenotype  is  a  result  of  or  response  to  castration. 

Milestone  1:  NEK6  can  confer  castrate  resistance  to  established  tumors,  kinase  activity  is  required,  and  continued 
expression  is  required.  The  NEK6-mediated  castration-resistant  tumors  are  AR  negative  and  squamous  in 
histology. 

Aim  2:  Assessing  signaling  pathways  involved  in  NEK6  and  NLK-mediated  castrate  resistance  and  assessing  their  clinical 
relevance. 

Milestone  2:  Determine  whether  AR  and  STAT3  are  necessary  to  confer  NEK6  and  NLK-mediated  castrate  resistance 

Milestone  3:  Determine  impact  ofNEK6  and  NLK  on  AR  and  STAT3  compared  to  other  hits  in  the  screen  in  vitro  and  in 
vivo 


Given  that  persistence  of  AR  signaling  has  been  demonstrated  to  be  an  important  mechanism  of  castration 
resistance,  we  sought  to  determine  if  NEK6  influences  AR  signaling.  As  noted  previously,  the  NEK6-mediated  castration 
resistant  tumors  lack  expression  of  the  AR,  making  it  unlikely  that  AR  is  responsible  for  this  phenotype.  NEK6 
overexpression  does  not  lead  to  an  increase  in  activity  of  an  AR  reporter  based  on  the  PSA  enhancer  in  FNCaP  cells 
(Supplemental  Figure  2A),  and  inducible  overexpression  of  NEK6  in  FHSR-AR  cells  (Supplemental  Figure  2B)  does  not 
lead  to  increased  expression  of  the  AR  targets  PSA  and  TMPRSS2  (Supplemental  Figure  2C).  Gene  set  enrichment 
analysis  (Subramanian  et  al.,  2005)  of  an  mRNA  expression  signature  of  NEK6  activity  generated  from  cells  inducibly 
expressing  wild-type  vs.  a  kinase  dead  NEK6  form  of  NEK6  (with  lysine  resides  at  positions  74  and  75  mutated  to 
methionine  -  O’Regan  and  Fry,  2009)  in  vitro  demonstrates  that  gene  expression  changes  associated  with  NEK6  kinase 
activity  do  not  correlate  positively  or  negatively  (Supplemental  Figure  2D)  with  three  previously  published  signatures  of 
AR  activity  (Hieronymus  et  al.,  2006;  Mendiratta  et  al.,  2009;  Sharma  et  al,  2013). 

It  has  previously  been  demonstrated  that  xenograft  tumors  with  squamous  histology  can  also  be  generated  by 
transformation  of  patient-derived  prostate  basal  cells  (Stoyanova  et  al.,  2013)  with  these  tumors  marked  by  active  |3- 
catenin;  however  there  is  no  evidence  for  nuclear  (3-catenin  staining  in  the  NEK6-mediated  castration  resistant  tumors 
generated  here  (Supplemental  Figure  3A).  Thus,  we  sought  to  determine  whether  other  canonical  oncogenic  signaling 
pathways  previously  reported  to  be  activated  by  NEK6  may  be  playing  a  role.  NEK6  purified  from  rat  liver  was  found  to 
be  the  major  protein  kinase  that  is  active  on  the  p70  S6  kinase  (RPS6KB1)  hydrophobic  regulatory  site,  Thr412  (Belham 
et  al.,  2001).  Subsequently,  it  was  demonstrated  that  NEK6  could  phosphorylate  hydrophobic  motifs  of  RPS6KB1  as  well 
as  SGKl  in  vitro  (Fizcano  et  al.,  2002).  In  addition,  STAT3  phosphorylation  at  serine  727  has  been  implicated  in  the 

6 


transformation  activity  of  NEK6  previously  (Jeon  et  aL,  2010).  However,  inducible  expression  of  NEK6  in  LHSR-AR 
cells  does  not  increase  phosphorylation  on  RPS6KB1  at  Thr412  or  STAT3  at  Ser727  with  or  without  growth  factor 
stimulation  as  compared  to  cells  expressing  kinase-dead  NEK6  or  GFP  (Supplemental  Figure  3B).  We  were  not  able  to 
find  a  quality  antibody  recognizing  Ser422  of  SGKl,  but  we  were  unable  to  demonstrate  phosphorylation  of  SGKl  at 
Ser422  by  active  recombinant  NEK6  in  an  in  vitro  kinase  assay  (data  not  shown). 

Because  NEK6  has  previously  been  described  to  play  a  role  in  the  G2/M  transition  (Yin  et  aL,  2003)  with  its 
proliferative  effects  in  hepatocellular  carcinoma  reported  to  be  mediated  through  modulation  of  cyclin  B  (Zhang  et  aL, 
2014),  we  assessed  whether  NEK6  overexpression  can  influence  cell  cycle  progression  in  our  system.  The  baseline  cell 
cycle  profile  is  identical  for  cells  with  inducible  expression  of  NEK6  with  and  without  doxycycline  in  normally  cycling 
cells,  cells  starved  in  growth-factor  free  media,  and  in  cells  released  into  growth-factor  containing  media  (Supplemental 
Figure  4A).  NEK6  does  not  increase  proliferation  rate  of  FHSR-AR  cells  in  vitro  compared  to  a  lacZ  control 
(Supplemental  Figure  4B).  Thus,  the  functional  role  of  NEK6  overexpression  in  these  cells  does  not  appear  to  be  related 
to  cell  cycle  progression.  Another  activity  of  NEK6  that  has  been  implicated  in  its  oncogenic  activity  is  inhibition  of  p53- 
mediated  senescence  (Jee  et  al.,  2010).  FHSR-AR  cells  express  large-T  antigen,  so  the  p53  pathway  should  be  inactive  in 
these  cells.  To  confirm  this  in  our  multiply  infected  FHSR-AR  cells,  they  were  exposed  to  etoposide  at  10  pM  for  18 
hours  or  50  pM  for  4  hours:  under  neither  condition  is  p21  expression  induced,  while  the  level  of  cleaved  PARP  is  not 
altered  by  NEK6  expression  (Supplemental  Figure  4C).  Thus  NEK6  does  not  act  by  modulating  the  p53  pathway,  and 
furthermore  the  pro-survival  effect  of  NEK6  appears  to  be  specific  to  the  in  vivo  context  of  castration  resistance. 

Milestones  2  and  3:  NEK6  does  not  activate  AR  signaling  and  does  not  lead  to  detectable  STAT3  phosphorylation 
in  LHSR-AR  cells.  In  addition,  there  is  no  evidence  that  NEK6  activates  beta-catenin  or  RPS6KB1  signaling, 
promotes  cell  cycle  progression,  or  antagonizes  p53  signaling  in  our  model. 

Milestone  4:  Generation  of  mRNA  and phosphoproteomic  signatures  corresponding  to  androgen  independence  conferred 
by  our  hits,  and  comparison  to  androgen  signaling  and  existing  databases 

Identification  of  FOXJ2,  NCOA5,  CSNKIAI  and  YESl  as  NEK6  substrates 

To  gain  a  more  comprehensive  understanding  of  the  immediate  signaling  mediated  by  NEK6  expression  and 
discover  novel  in  vivo  substrates,  we  performed  a  phosphoproteomic  analysis.  Constructs  for  the  expression  of  wild-type 
and  kinase-dead  NEK6  under  the  control  of  a  doxycycline-inducible  promoter  were  introduced  into  LHSR-AR  cells,  and 
we  created  phosphoproteomic  profiles  of  cells  expressing  doxycycline  inducible  versions  of:  (1)  wild- type  NEK6  in  the 
presence  of  doxycycline,  (2)  kinase-dead  NEK6  in  the  presence  of  doxycycline,  and  (3)  wild  type  NEK6  without 
doxycycline.  Cells  were  cultured  in  “light”,  “medium”,  and  “heavy”  media  for  Stable  Isotope  Labeling  by  Amino  acids  in 
Cell  culture  (SILAC)  for  7  days  in  2  different  permutations;  they  were  growth  factor  starved  for  the  final  30  hours  and 
then  growth  factor  stimulated  for  5  minutes.  Lysates  were  then  subjected  to  SCX-IMAC  phosphorylation  analysis  and 
assayed  by  mass  spectrometry. 

A  total  of  9418  phosphopeptides  (8432  phosphoserine,  952  phosphothreonine,  34  phosphotyrosine)  from  3401 
proteins  were  detected  in  this  experiment.  59  phosphopeptides  from  50  proteins  were  increased  in  both  the  wild- type 
induced  vs.  un-induced  and  wild-type  vs.  kinase  dead  comparison  with  a  combined  q  value  of  <0.25  (Supplemental  Table 
2).  The  differentially  phosphorylated  proteins  include  a  large  number  of  transcriptional  regulators,  including  the 
forkhead-box  family  proteins  FOX03  and  FOXAl,  and  the  phosphopeptides  represent  a  variety  of  common 
phophorylation  motifs  as  shown  in  Table  2.  Among  these  motifs  are  pS/pT-P,  associated  with  MAPK/CDK/GSK3 
signaling,  and  R-X-R-X-X-pS/pT,  associated  with  AKT  and  RSK  (ribosomal  S6-kinase)  signaling.  We  noted  that  the 
frequency  of  these  motifs  among  the  59  differential  phosphopeptides  is  not  increased  as  compared  to  their  frequency 
among  all  detected  phosphopeptides  -  for  example,  pS/pT-P  motifs  were  found  in  26  of  the  59  differential 
phosphopeptides  vs.  4815  of  the  9434  total  phosphopeptides  (44%  vs.  51%,  p=NS).  This  suggests  that  NEK6  does  not 
globally  activate  MAPK/CDK/GSK3,  AKT/RSK,  CAMK  or  CK2  signaling. 

However,  a  canonical  NEK6  motif  (Vaz  Meirelles  et  al.,  2010)  with  leucine  at  the  -3  position  [L-X-X-pS/pT- 
F/W/Y/M/L/I/V/R/K]  is  more  frequently  represented  in  the  59  significantly  enriched  phosphopeptides  compared  to  all 
detected  phosphopeptides  (p=0.0092  Fisher’s  exact);  an  “acceptable”  NEK6  substrate  motif  [L/F/W/Y-X-X-pS/pT- 
F/W/Y/M/L/I/V/R/K]  is  even  more  significantly  enriched  (p=  1.95x10'^).  This  suggests  that  this  peptide  sequence  is  a  true 
description  of  the  NEK6  phosphorylation  motif,  and  that  the  8  proteins  with  phosphopeptides  of  this  form  detected  here 
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(FOXJ2,  HUWEl,  NCOA5,  KRT18,  TRA2B,  HNRNPM,  HNRNPA2B1,  ZNF326)  are  bona  fide  in  vivo  NEK6 
substrates. 


Table  2.  Summary  of  phosphoproteomic  data 

L-X-X-pS/pT- 

F/W/Y/M/L/I/ 

V/R/K 

(NEK6  general 
eonsensus) 

E/F/W/Y-X-X- 

pS/pT- 

F/W/Y/M/E/I/ 

V/R/K 

(NEK6 

acceptable) 

pS/pT-P 
(pro  line- 
directed 
kinases) 

r_X-R-X-X- 
pS/pT 
(AKTl  or 
RSK  family) 

R-X-X-pS/pT 

(CaMK 

consensus) 

pS/pT-D/E-X- 

D/E 

(CK2 

consensus) 

Others 

Proteins  with 
phosphopeptides 
meeting  signifieanee 
threshold 
(59  total 

phosphopeptides 
representing  50 
proteins) 

FOXJ2 

HUWEl 

NCOA5 

KRT18 

FOXJ2 

HUWEl 

NCOA5 

KRT18 

TRA2B 

(Y-X-X-pS-Y) 

HNRNPM 

(F-X-X-pS-F) 

HNRNPA2B1 

(F-X-X-pS-F) 

ZNF326 

(F-X-X-pS-Y) 

FOX03 

TRPSl 

LM07 

SATB2 

PAK6 

BCE6 

FOXAl 

KLF4 

ATXNl 

IRF2BP1 

RFX2 

HIVEP2 

KEF3 

Others  (xl3) 

SLC2A12 

FAM21C 

ATXNl 

SLC2A12 

FAM21C 

ATXNl 

ATM 

LM07 

DROSHA 

PEEKHA6 

RIPK3 

MYOF 

EXPH5 

MLET3 

ERCC5 

LIGl 

PBRMl 

16 

phosphopeptid 
es  representing 
14  proteins 

12  pS 

3pT 

1  pY 

Enriehment 
(eompared  to  all 
deteeted 

phosphopeptides)^ 

4/59  vs. 
131/9362 
(p=0.0092 
Fisher’s  exaet) 

8/59  vs. 
186/9362 
(p=  1.95x10'^) 

26/59  vs. 
4815/9434 
(p=NS) 

3/59  vs. 
434/9301 
(P=NS) 

11/59  vs. 
2239/9362 
(p=NS) 

3/59  vs. 
701/9417 
(P=NS) 

^  The  denominator  1 
phosphopeptides  ab 
detected  by  mass  sp 

for  the  total  number  of  detected  phosphopeptides  in  each  column  is  filtered  for  only  those 

fe  to  be  defined  as  representing  this  motif  (i.e.  only  those  phosphopeptides  where  the  -5  amino  acid  was 

lectrometry  were  assessed  for  representing  the  AKT/RSK  motif) 

To  confirm  whether  detection  of  these  phosphopeptides  in  cells  is  due  to  direct  phosphorylation  by  NEK6,  we 
assessed  the  activity  of  commercially  available  recombinant  GST-NEK6  on  a  subset  of  these  proteins  immunoprecipitated 
via  a  C-terminal  V5-tag  through  an  on-bead  kinase  assay.  We  initially  focused  on  the  transcription  factors  FOXJ2  and 
NCOA5  because  they  are  related  to  families  of  transcription  factors  already  described  to  be  important  in  prostate  cancer 
(Tao,  et  aL,  2014;  Qin  et  aL,  2014).  We  were  able  to  confirm  phosphorylation  of  FOXJ2  and  NCOA5,  with  no 
phosphorylation  of  a  lacZ  control  (Figure  5A).  Mutation  of  the  phosphorylation  sites  of  FOXJ2  and  NCOA5  identified 
via  mass  spectrometry,  Ser8  and  Ser96  respectively,  decreases  the  in  vitro  phosphorylation  of  these  substrates, 
demonstrating  the  specificity  of  the  kinase  activity  at  these  residues. 


To  expand  our  phosphoproteomic  studies  in  vivo  we  sought  to  elucidate  NEK6-interacting  proteins  and  potential 
substrates  in  the  tumors  themselves.  To  identify  in  vivo  interacting  proteins  using  mass  spectrometry,  we  analyzed  tumors 
that  formed  in  female  mice  from  LHSR-AR  cells  overexpressing  a  form  of  NEK6  with  a  C-terminal  V5  tag  or 
overexpressing  untagged  NEK6;  protein  lysates  from  these  tumors  were  subjected  to  immunoprecipitation  with  anti-V5 
magnetic  beads.  Proteins  enriched  in  the  immunoprecipitate  from  NEK6-V5  expressing  tumors  compared  to  tumors 
expressing  NEK6  without  a  V5-tag  (from  which  the  immunoprecipitate  would  consist  only  of  proteins  bound  non- 
specifically  to  the  beads)  were  considered  as  candidate  in  vivo  interacting  proteins.  A  total  of  258  proteins  were  identified 
that  were  enriched  with  a  log2  ratio  of  greater  than  1.25  in  the  experimental  sample  vs.  the  control  (Supplemental  Figure 
3).  Gene  Ontology  (GO)  terms  associated  with  this  set  of  proteins  as  determined  through  the  Gene  Ontology  for 
Functional  Analysis  (GOFF A)  software  (Sun  et  aL,  2006)  are  listed  in  Table  3. 


From  this  list,  we  focused  on  the  interacting  protein  kinases,  as  direct  substrates  of  NEK6  that  are  protein  kinases 
would  likely  mediate  further  signaling  downstream  from  NEK6  overexpression.  We  obtained  ORFs  for  nine  of  the 
interacting  proteins  identified  in  this  analysis  and  then  performed  an  in  vitro  kinase  assay  with  NEK6.  As  seen  in  Figure 
5A,  NEK6  potently  phosphorylates  the  kinases  CSNKIAI  and  YESl  in  vitro. 


Table  3.  NEK6-interacting  proteins  from  tumors  categorized  1 

by  GO  term 

Cytoskeleton 

Nucleoside 

Calcium 

Signal  transduction 

Protein 

Nucleotid 

Misc. 

triphosphatase  activity 

ion 

kinase 

e 

binding 

metabolic 

process 

ACTCl 

MYH9 

ACTCl 

MYOID 

ACTN4 

ALDHIA 

NPMl 

CAD 

ACTCl 

Ribosome 

ACTN4 

MYL12B 

ATAD3B 

MYOIE 

AIFIL 

1 

PHB 

CAMK2 

ALDHIA 

GNB2L1 

ACTRIA 

MYL6 

ATL3 

MY05A 

ANXAl 

ANXAl 

PPP2R2A 

D 

1 

RPLIOA 

ACTR2 

MYL9 

ATPlAl 

MY05C 

ANXA2 

(CALMl 

PRDX4 

CPNE3 

ATL3 

RPL14 

ACTR3 

MY018 

ATP5C1 

MY06 

(CALMl) 

) 

PSMCl 

CSNKIA 

ATPlAl 

RPL18 

AIFIL 

A 

ATP50 

PSMCl 

CALML3 

CAMK2 

PSMC5 

1 

ATP1A4 

RPL18A 

ALDOC 

MYOIB 

DCTNl 

PSMC5 

GSN 

D 

PSMC6 

CSNK2B 

ATP5C1 

RPL27 

ANXAl 

MYOID 

DDX3X 

PSMC6 

MYL12B 

CDC37 

PSMDl 

IKBKB 

ATP50 

RPL4 

ARPCIA 

MYOIE 

EEFlAl 

RABllB 

MYL6 

CLTA 

PSMD14 

MAP2K2 

CAD 

RPL7A 

ARPC2 

MY05A 

EEF2 

RAB14 

MYL9 

CNGBl 

RABllB 

MAPKl 

CNP 

RPS15A 

ARPC4 

MY05C 

EIF4A1 

RAB35 

PLS3 

CSNKIA 

RAB14 

NEK6 

CTPS 

RPS25 

ARPC5 

MY06 

EIF4A3 

RAB5A 

SIOOAIO 

1 

RABIB 

NEK7 

DOCK7 

RPS27 

ARPC5L 

NPMl 

GNA13 

RAB7A 

SlOOAll 

CSNK2B 

RAB35 

NEK9 

FLNA 

RPS3 

(CALMl 

PDCD6I 

GNBl 

RAN 

S100A2 

CTNNDl 

RAB5A 

NME2 

GNA13 

RPS3A 

) 

P 

GNB2 

RHOC 

S100A4 

DDBl 

RAB7A 

PRKDC 

GNBl 

RPS4X 

CAPZAl 

PDLIMl 

HSP90A 

RRAS2 

S100A8 

DOCK7 

RAN 

RPS6KA 

GNB2 

RPS5 

CAPZB 

PDLIM5 

A1 

RUVBLl 

S100A9 

DPYSL2 

RHOC 

I 

GNB2L1 

UBA52 

CCT3 

POTEKP 

HSPA8 

RUVBL2 

SLC25A1 

FLNA 

RPS6KA 

RPS6KA 

HSP90A 

CCT4 

RAI14 

KIF21B 

SETX 

3 

FLNB 

1 

3 

A1 

Proteaso 

CCT5 

RCC2 

KIF5B 

TUBA  1C 

SSR4 

GMFB 

RPS6KA 

TWFl 

IQGAPl 

mal 

CDK5RA 

RUVBLl 

KRAS 

TUBB 

TPM4 

GNA13 

3 

YESl 

MYH9 

complex 

P2 

S100A8 

MCM5 

TUBB2A 

GNBl 

RRAS2 

NAMPT 

PSMCl 

CEP350 

S100A9 

MYHIO 

TUBB2C 

GNB2 

SIOOAIO 

NME2 

PSMC5 

CNN2 

SEPT2 

MYH14 

TUBB6 

GNB2L1 

SlOOAll 

NPPC 

PSMC6 

CORO  1 A 

SPRRIB 

MYH9 

TUFM 

GSN 

S100A4 

PSMCl 

PSMDl 

CORO  IB 

SPTBN2 

MYL6 

VCP 

GSTPl 

S100A9 

PSMC5 

PSMDl  4 

CSNKIA 

SVIL 

MY018A 

YMEILI 

HNRNP 

SFN 

PSMC6 

VCP 

1 

TCPl 

MYOIB 

K 

TOLLIP 

RABllB 

CTTN 

TMODl 

HSP90A 

UACA 

RAB14 

Unfolded 

DBT 

TMOD3 

A1 

UBA52 

RAB35 

protein 

DCTNl 

TPMl 

IKBKB 

VAPA 

RAB5A 

binding 

DSTN 

TPM4 

IQGAPl 

VCP 

RAB7A 

CALR 

FLNA 

TUBA  1C 

KRAS 

YWHAE 

RAN 

CCT3 

FLNB 

TUBB 

KRT18 

YWHAG 

REX02 

CDC37 

GNB2L1 

TUBB2A 

MAP2K2 

YWHAH 

RHOC 

DNAJA2 

GSN 

TUBB2C 

MAPKl 

YWHAZ 

RRAS2 

DNAJA4 

IQGAPl 

TUBB6 

MYH9 

RUVBL2 

HSP90A 

KIF21B 

TWFl 

MYOIE 

SLC25A1 

A1 

KIF5B 

UACA 

MY06 

3 

HSPA2 

KRT18 

YWHAE 

NAMPT 

SULT2B 

HSPA9 

KRT8 

YWHAZ 

NEK6 

1 

NAP1L4 

LIMAl 

ZNF185 

TPIl 

PPIA 

MYHIO 

TPMl 

PPIB 

MYH14 

VCP 

XDH 

While  the  two  approaches  of  in  vitro  phosphoproteomic  analysis  using  SILAC  and  in  vivo  interactomic  analysis 
can  identify  a  large  number  of  putative  substrates,  neither  approach  is  comprehensive  due  to  limitations  of  the  detection 
platforms.  Thus,  we  sought  to  validate  and  determine  the  biological  consequences  of  NEK6  phosphorylation  of  the  top 
candidates  identified  through  both  approaches. 
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Mapping  NEK6  phosphorylation  sites  on  FOXJ2,  NCOA5,  CSNKIAI  and  YESl 


To  comprehensively  map  the  phosphorylation  sites  of  NEK6  on  FOXJ2,  NCOA5,  CSNKIAI  and  YESl,  all 
serine  and  threonine  residues  in  these  proteins  in  “acceptable”  NEK6  phosphorylation  motifs  as  described  above 
[L/F/W/Y-X-X-pS/pT-F/W/Y/M/L/I/V/R/K]  were  identified;  when  the  relevant  serines  and  threonines  are  mutated  to 
aspartic  acid,  the  resulting  mutant  forms  of  these  proteins  are  phosphorylated  significantly  less  than  the  wild  type 
counterpart  (compare  lanes  2  and  1  in  right  panels  of  Figures  5B  and  C,  and  Figures  6B  and  C).  When  mutating  all  of 
these  potential  phosphorylation  sites  but  maintaining  one  of  the  serines/threonines  as  wild  type,  it  is  apparent  that  NEK6 
can  phosphorylate  FOXJ2  at  Thr23  and  Ser254  (Figure  5B,  lanes  3  and  5)  in  addition  to  Ser8;  NCOA5  at  Ser21  and 
SerlSl  (Figure  5C,  lanes  3  and  4)  in  addition  to  Ser96;  CSNKIAI  at  Ser96,  Thrl29,  and  Thr287  (Figure  6B,  lanes  4-6) 
and  YESl  at  Thr295  and  Ser524  (Figure  6C,  lanes  5-6). 

Gene  Expression  changes  mediated  by  NEK6 

To  further  elucidate  downstream  signaling  mediated  by  NEK6  overexpression,  we  sought  to  determine  gene 
expression  changes  in  tumors  with  NEK6  overexpression  in  castrate  conditions.  We  established  tumors  at  3  subcutaneous 
sites  of  male  mice  implanted  with  a  testosterone  pellet  from  cells  expressing  NEK6  under  a  doxycycline-inducible 
promoter.  After  tumors  were  formed,  doxycycline  was  withdrawn  from  half  the  mice  7  days  prior  to  planned  tumor 
harvest.  Mice  were  then  castrated  and  testosterone  pellets  removed;  they  were  then  sacrificed  and  tumors  were  harvested 
at  days  0,  2  and  5  after  castration.  As  seen  in  Figure  7 A,  withdrawal  of  doxycycline  reduces  NEK6  levels  to  near  the 
levels  in  tumors  without  exogenous  NEK6  expression.  Messenger  RNA  was  then  isolated  from  tumors  from  the  day  2 
and  day  5  time  points  and  subjected  to  gene  expression  analysis  using  RNASeq. 

Gene  expression  changes  with  NEK6  overexpression  in  tumors  at  day  5  after  castration  demonstrate  some  overlap 
with  upregulated  and  downregulated  genes  in  the  comparison  of  NEK6  wild-type  vs.  kinase-dead  overexpression  in 
LHSR-AR  cells  in  vitro  in  the  experiment  described  in  Supplemental  Figure  ID  (Figure  7B).  Given  that  the  gene 
expression  was  measured  in  different  contexts  (tumors  vs.  culture)  and  using  different  assays  (RNASeq  vs.  Affymetrix 
microarray)  in  both  comparisons,  the  genes  in  this  overlap  are  high  confidence  NEK6  downstream  transcriptional  targets. 
Interestingly,  the  upregulated  genes  include  prostate  stem  cell  antigen  (PSCA),  which  is  a  marker  of  late  intermediate 
prostate  epithelial  cells  (Tran  et  aL,  2003)  that  has  been  reported  to  be  increased  in  prostate  cancer  metastasis  (Lam,  et  al. 
2005);  and  keratin  13  (KRT13),  which  is  a  marker  of  squamous  differentiation  in  the  pathology  literature  (van  Dorst  et  al, 
1998). 


To  further  characterize  these  gene  expression  changes.  Gene  Set  Enrichment  Analysis  was  performed  on  the 
signature  derived  from  tumors  overexpressing  NEK6  at  day  5  after  castration  compared  to  controls  (Figure  7C).  The  top 
two  curated  gene  sets  (C2)  from  the  molecular  signatures  database  (MSigDB)  correlated  with  the  NEK6  signature  are  1) 
genes  down-regulated  in  primary  B  lymphocytes  within  60-180  min  after  activation  of  LMPl,  an  oncogene  encoded  by 
Epstein-Barr  virus  (Dirmeier  et  al,  2005)  and  2)  genes  representing  the  epithelial  differentiation  module  in  sputum  during 
asthma  exacerbations  (Bosco,  et  aL,  2010).  GOFFA  reveals  the  most  statistically  significant  GO  biological  process 
associated  with  NEK6  upregulated  genes  to  be  “response  to  biotic  stimulus”  (p=0.0005);  other  relevant  terms  include 
“type  I  interferon  signaling  pathway”  (p=0.0023),  “apoptotic  signaling  pathway”  (p=0.0056),  “cytoskeleton  organization” 
(p=0.0084)  and  “epithelial  cell  differentiation”  (p=0.0183).  More  broadly,  68  NEK6  upregulated  genes  are  categorized 
under  “developmental  process,”  and  33  genes  are  categorized  under  “immune  system  process,”  so  we  used  these  GO 
terms  as  a  framework  for  subsequent  assessments  of  gene  expression  changes  in  tumors. 

We  hypothesized  that  NEK6  could  mediate  resistance  to  castration  in  tumors  by  maintaining  survival  signaling 
that  is  lost  when  AR  is  no  longer  activated  by  circulating  androgens.  To  understand  components  of  this  survival 
signaling,  we  compared  gene  expression  of  control  tumors  (i.e.  without  NEK6  overexpression)  at  the  day  2  and  day  5  time 
points.  Using  gene  set  enrichment  analysis  (GSEA)  we  find  that  the  top  gene  set  and  five  of  the  top  10  gene  sets 
associated  with  gene  expression  lost  following  castration  reflect  interferon  signaling  (Table  4,  Figure  8A  left).  This 
suggests  that  in  our  in  vivo  model,  much  of  the  survival  signaling  mediated  by  circulating  androgens  involves  expression 
of  genes  associated  with  these  pathways.  Interferon-related  gene  sets  have  previously  been  described  as  being 
biologically  important  in  prostate  cancer,  particularly  with  regards  to  differences  in  signaling  between  cancers  in  patients 
with  African  and  European  ancestry  (Wallace  et  al.,  2008  -  Figure  8A,  right)  and  in  treatment  resistance  (Weichselbaum, 
et  al.,  2008;  Cheon  et  al.,  2013). 
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Table  4.  Gene  sets  enriched  in  ranked  list  of  genes  with  decreased  expression  after  castration  (Control  Day  5  vs.  Day  2) 


GS 

SIX 

ES 

NES 

NOM 

p=val 

FDR 

q-val 

FWE 

Rp- 

val 

1 

REACTOMEJNTERFERON_ALPHA_BETA_SIGNALING 

46 

0.69 

3.2 

0 

0 

0 

2 

LIANG_SILENCED_BY_METHYLATION_2 

44 

0.64 

3.04 

0 

0 

0.001 

3 

BOWIE_RESPONSE_TO_TAMOXIFEN 

17 

0.82 

3 

0 

0.001 

0.003 

4 

BROWNEJNTERFERON_RESPONSIVE_GENES 

61 

0.6 

2.98 

0 

0.001 

0.005 

5 

ZHANGJNTERFERON_RESPONSE 

20 

0.78 

2.98 

0 

0.001 

0.005 

6 

BOWIE_RESPONSE_TO_EXTRACELLULAR_MATRIX 

16 

0.8 

2.97 

0 

0.001 

0.009 

7 

MOSERLEJFNA_RESPONSE 

30 

0.68 

2.93 

0 

0.002 

0.013 

8 

SANA_TNF_SIGNALING_UP 

71 

0.56 

2.92 

0 

0.002 

0.016 

9 

LI  U_VAV3_PROSTATE_CARCI  NOG  EN  ESIS_U  P 

59 

0.57 

2.88 

0 

0.003 

0.023 

10 

DERJFN_ALPHA_RESPONSE_UP 

69 

0.55 

2.88 

0 

0.002 

0.023 

Gene  Expression  changes  mediated  by  NEK6  substrates 

In  order  to  determine  the  role  of  the  identified  substrates  in  mediating  signaling  downstream  of  NEK6,  we  assess 
signaling  in  tumors  mediated  by  phosphomimic  forms  of  these  substrates  with  serine  to  aspartic  acid  substitutions  of  the 
phosphorylation  sites  identified  in  vitro.  These  phosphomimic  forms  were  not  sufficient  either  alone  or  in  combination  to 
generate  castration-resistant  tumors  (data  not  shown).  In  order  to  interrogate  downstream  signaling  mediated  by  these 
substrates  in  castration-resistant  tumors,  we  expressed  them  under  a  doxycycline-inducible  promoter  in  tumors  formed  in 
male  mice,  which  were  harvested  for  RNASeq  five  days  after  castration  from  mice  with  continued  doxycycline  treatment 
or  7  day  doxycycline  withdrawal  as  for  NEK6  previously.  GSEA  reveals  that  the  gene  expression  changes  mediated  by 
phosphomimic  forms  of  CKla,  FOXJ2,  and  YESl,  but  not  NCOA5,  are  statistically  significantly  correlated  with  the 
NEK6  signature.  These  correlations  are  not  due  to  technical  artifact  or  systemic  bias  in  the  tumor  samples  with  continued 
expression  of  the  transgene  (+dox)  vs.  discontinuation  of  expression  (-dox)  because  1)  the  genes  upregulated  by 
expression  of  phosphomimic  forms  of  FOXJ2,  CKla  and  YESl  are  largely  non-overlapping  2)  a  transgene  assayed  in  the 
same  system,  NCOA5,  does  not  demonstrate  statistically  significant  overlap  and  3)  when  performing  comparative  marker 
selection  for  all  +dox  samples  vs.  the  -dox  samples,  no  gene  is  upregulated  in  the  comparison  with  a  signal  to  noise  ratio 
greater  than  1 . 

The  phosphomimic  substrate  with  the  highest  degree  of  correlation  with  the  NEK6  signature  is  CKla  (Figure  8C). 
The  set  of  genes  that  are  upregulated  by  both  CKla(3S^D)  and  NEK6  is  highly  enriched  for  genes  involved  in 
developmental  processes  (the  top  GO  biological  process  with  >3  genes  in  the  overlap  is  “epithelial  cell  differentiation” 
p=0.0012).  Included  in  this  overlap  is  KRT13,  a  marker  of  squamous  differentiation,  suggesting  that  CKla  is  involved  in 
the  transdifferentiation  of  NEK6-mediated  castration  resistant  tumors  to  a  squamous  phenotype.  In  addition,  RNASeq 
suggests  that  FOXJ2  is  the  relevant  transcription  factor  for  some  of  the  high  confidence  NEK6  transcriptional  targets 
(TPPP3,  PSCA,  PLAC8)  from  Figure  7B,  and  that  YESl  is  also  involved  in  aspects  of  the  differentiation  phenotype 
(Supplemental  Figure  5). 

As  mentioned  previously,  transduction  of  LHSR-AR  cells  with  a  combination  of  phosphomimic  forms  of  the 
identified  substrates  is  not  sufficient  to  lead  to  androgen-independent  tumor  formation.  Interestingly,  none  of  these 
phosphomimics  leads  to  much  upregulation  of  “immune  response”  genes,  and  we  hypothesize  that  this  signaling  may  be 
important  for  the  castration  resistance  phenotype,  particularly  since  the  genes  downregulated  by  castration  in  this  model 
are  highly  correlated  with  interferon/TNF  signaling.  The  NEK6  interactome  as  shown  in  Table  3  includes  many  genes 
putatively  involved  in  immune  signaling  including  IKBKB,  DDX3X,  ILF2,  PRDX4,  RHOC,  S100A8,  S100A9,  and 
TOLLIP;  however,  none  of  these  proteins  act  as  direct  NEK6  substrates  in  vitro  in  comparison  to  the  potent  substrates 
CKla  and  YESl  (Supplemental  Figure  6).  In  addition,  NEK6  has  previously  been  reported  to  interact  with  RELB  (Vaz 
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Meirelles  et  aL,  2010)  and  increase  the  activity  of  an  NFkB  promoter  (Matsuda  et  aL,  2003),  and  it  has  also  been 
previously  reported  to  phosphorylate  STATS  (Jeon  et  al.,  2010);  however,  RELA,  RELB,  STATl,  STATS,  and  STATS 
are  also  not  phosphorylated  to  an  appreciable  extent  by  NEK6  in  our  in  vitro  kinase  assay.  Thus,  the  impact  of  NEK6  on 
immune  signaling  may  be  indirect  or  the  relevant  substrate(s)  have  yet  to  be  identified. 

Milestone  4:  Phosphoproteomic  studies  reveal  novel  substrates  FOXJ2,  HUWEl,  NCOA5,  KRT18,  TRA2B, 
HNRNPM,  HNRNPA2B1,  ZNF326.  Proteomic  studies  identify  a  number  of  potential  interacting  proteins  in 
tumors;  of  these  CSNKIAI  and  YESl  are  also  in  vitro  substrates.  Gene  expression  studies  reveal  NEK6  likely 
mediates  castration  resistance  through  alteration  of  cell  differentiation  state  and  maintenance  of  immune-related 
signaling  that  is  normally  lost  with  castration. 

Milestone  5:  Determine  evidence  for  NEK6  and  NLK  conferring  castrate  resistance  in  clinical  samples  and  identify 
subsets  of  patients  displaying  genomic  changes  and  expression  signatures  corresponding  to  resistance  mediated  through 
the  activity  of  a  particular  gene/pathway. 

NEK6  is  overexpressed  in  several  malignant  tissues  and  cell  lines,  and  has  been  previously  been  implicated  in  cell 
transformation  (Jeon  et  ai,  2010;  Nassirpour  et  al.,  2010;  Wang  et  aL,  2014).  Cancer  types  previously  described  to  have 
high  frequency  of  NEK6  overexpression  include  cancers  of  the  liver  (Chen  et  ai,  2006;  Jeon,  et  aL,  2010;  Cao,  et  ai, 
2012),  stomach  (Takeno,  et  al.,  2008;  Nassirpour  et  al.,  2010),  breast  (Jeon  et  aL,  2010;  Nassirpour  et  al.,  2010),  uterus, 
colon/rectum  and  ovary  (Nassirpour  et  al.,  2010).  In  human  prostate  cancer,  the  NEK6  locus  at  chromosome  9q33.3  is 
present  in  a  region  of  low  level  copy  number  gain  (Taylor  et  ai,  2010;  Huang  et  ai,  2012;  Grasso  et  al.,  2012)  without  a 
known  validated  prostate  cancer  oncogene.  Genomic  characterization  of  human  prostate  cancers  through  the  Prostate 
TCGA  and  SU2C-PCF  Dream  Team  efforts  demonstrates  amplification  or  mRNA  overexpression  (z-score>2)  of  NEK6  in 
7%  of  primary  tumors  and  6%  of  metastases  (www.cbioportal.org/prostate-portal/). 

We  performed  immunohistochemistry  analysis  of  primary  prostate  cancers  compiled  in  seven  tumor  microarrays 
(TMA)  generated  by  the  Gelb  Center  for  Translation  Research  at  Dana-Farber  Cancer  Institute  using  an  antibody  we 
optimized  for  IHC  (Supplemental  Figure  2A  and  B).  NEK6  expression  was  quantified  by  spectral  imaging.  We  failed  to 
find  detectable  levels  of  NEK6  protein  in  most  prostate  tissues,  but  found  that  NEK6  is  aberrantly  expressed  in  -^16%  of 
tumors  (Supplemental  Figure  2C).  NEK6  expression  was  greater  in  tumor  samples  than  in  benign-appearing  glandular 
epithelium  from  these  same  radical  prostatectomy  specimens  (one-tailed  p=0.0012  by  Student’s  T-test),  suggesting  that 
NEK6  expression  is  associated  with  tumor  progression  (Figure  3A).  We  are  in  the  process  of  correlating  NEK6  status  with 
grade  and  outcomes  based  on  the  schema  in  Figure  3B. 

It  remains  unclear  at  this  point  whether  de  novo  acquisition  of  NEK6  amplification  or  overexpression  is  a 
mechanism  of  development  of  castration  resistance  in  advanced  prostate  cancer.  The  strongest  evidence  for  this  would  be 
detection  of  markers  of  overactive  NEK6  signaling  in  metastatic  CRPC  in  comparison  to  therapy-naive  samples,  which 
per  our  model  would  be  most  likely  to  be  detected  in  prostate  cancers  that  have  become  AR-independent.  A  recent  study 
by  Grasso  et  al  (2012)  genetically  profiled  35  cases  of  CRPC  in  comparison  to  59  cases  of  localized  PrCa.  The  primary 
difference  in  the  genomic  landscape  in  these  two  states  was  massive  amplification  of  the  AR  in  CRPC  in  comparison  to 
the  primary  cases;  there  was  no  increase  in  amplitude  of  the  region  of  copy  number  gain  on  9q33.3  where  NEK6  is 
located.  There  was  no  evidence  for  enrichment  in  genetic  markers  of  any  other  signaling  pathways  in  CRPC  cases, 
though  this  is  likely  due  to  the  fact  that  these  cases  were  collected  at  a  time  when  many  of  the  current  novel  hormonal 
therapies  had  not  yet  been  in  widespread  use.  Efforts  to  genetically  characterize  larger  numbers  of  CRPC  cases,  including 
those  from  patients  who  have  progressed  on  novel  hormonal  therapies,  are  currently  underway  at  DFCI 
(http://www.aacr.org/home/public-media/stand-up-to-cancer/su2c-dream-teams/su2c-pcf-dream-team-precision-therapy- 
of-advanced-prostate-cancer.aspx). 

We  analyzed  the  levels  of  NEK6  expression  in  several  patient-derived  prostate  cell  lines  by  immunoblotting 
(Figure  4 A)  and  found  that  most  prostate  cancer  cell  lines  expressed  higher  levels  of  NEK6  than  immortalized  (RWPE, 
PrEC-LH)  and  transformed  (LHSR-AR)  prostate  epithelial  cells.  Interestingly,  NEK6  levels  were  relatively  high  in  VCaP 
and  LNCaP  cells,  suggesting  that  high  expression  of  NEK6  was  not  sufficient  to  overcome  the  in  vitro  androgen 
dependence  of  these  cells;  however,  it  is  important  to  note  that  these  cells  were  derived  from  patients  with  castration 
resistant  disease  in  vivo.  To  test  whether  NEK6  was  essential  for  in  vitro  proliferation  of  androgen-independent  AR-low 
cell  lines  (reminiscent  of  the  AR-low  xenograft  tumors),  we  introduced  doxycycline-inducible  shRNAs  targeting  NEK6  to 
CL-1,  PC-3  and  DU145  cells  (Figure  4B).  Despite  near  complete  suppression  of  NEK6  expression,  there  was  no  effect  on 
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proliferation  of  the  cell  lines  except  a  modest  suppression  of  DU  145  proliferation  with  shRNA#2.  These  observations  are 
in  consonance  with  a  prior  report  suggesting  that  NEK6  is  not  essential  for  proliferation  of  many  mammalian  cell  lines 
(Nassirpour,  et  al.,  2010). 


Milestone  5:  NEK6  is  aberrantly  expressed  in  a  subset  of  primary  human  prostate  cancer.  Efforts  to  characterize 
metastatic  biopsies  from  patients  with  CRPC  for  sequencing  and  assessing  gene  expression  are  underway. 

Task  9:  Genome  wide  screen  of  genes  conferring  androgen-independence  to  androgen-dependent  cell  line 

Milestone  6:  Identify  a  set  of  genes  leading  to  androgen-independent  prostate  cancer  cell  growth  in  vitro  and 
in  vivo. 


To  ensure  adequate  representation  of  each  gene  in  the  kinome  library  in  our  original  in  vivo  functional 
genomic  screen,  we  implanted  mice  with  LHSR-AR  cells  transduced  with  9-10  ORF  constructs  in  each  pool.  It 
would  be  impractical  to  expand  this  screen  to  a  genome-wide  library  because  of  the  number  of  mice  and  amount 
of  material,  time  and  effort  required.  In  order  to  obtain  a  more  manageable  list  of  genes  to  query  in  our  in  vivo 
assay,  we  first  aim  to  identify  genes  that  can  confer  androgen-independent  proliferation  to  the  androgen- 
dependent  LNCaP  prostate  cancer  cell  line  in  vitro.  This  experiment  will  be  performed  as  a  pooled  screen  using 
a  barcoded  lentivirally-delivered  library  encompassing  ~20,000  ORFs  available  through  the  Broad  Institute 
Genetic  Perturbation  Platform  using  the  schema  in  Figure  9.  The  cells  will  be  infected  at  a  target  infection  rate 
of  30%  to  achieve  a  multiplicity  of  infection  (MOI)  of  1  and  target  representation  of  1000  cells  infected  with 
each  ORF.  The  representation  of  ORFs  at  each  time  point  will  be  assessed  by  isolating  genomic  DNA, 
performing  targeted  sequencing  of  the  barcodes,  and  deconvoluting  enriched  ORFs  as  represented  by  their 
corrensponding  barcodes.  ORF  representation  will  be  assessed  at  an  early  time  point  and  after  7  weeks  of 
selection,  either  in  media  with  androgen-poor  charcoal  stripped  serum  (CSS)  or  media+CSS+2.5  uM 
enzalutamide,  both  compared  with  cells  grown  in  androgen-containing  media  with  fetal  bovine  serum  (FBS)  as 
a  control.  Three  replicates  will  be  obtained  for  each  condition,  and  candidates  identified  through  the  screen  will 
be  validated  in  low  throughput,  and  also  assessed  in  vivo  for  conferring  androgen  independent  xenograft  tumor 
formation  to  LHSR-AR  cells  and  enzalutamide  resistance  to  xenografts  derived  from  LNCaP  cells 
overexpressing  the  androgen  receptor  (LNCaP-AR  cells). 

Conditions  for  achieving  infection  rate  of  30%  have  been  defined,  and  we  will  initiate  the  screen  over 
the  next  few  weeks. 

KEY  RESEARCH  ACCOMPLISHMENTS 

•  Identification  of  NEK6  as  a  novel  mediator  of  castration  resistance  in  an  in  vivo  forward  functional  genomic 
screen. 

•  Demonstration  that  kinase  activity  and  continued  expression  of  NEK6  is  required  for  maintenance  of  castration 
resistance,  suggesting  its  suitability  as  a  therapeutic  target. 

•  Discovery  of  novel  NEK6  substrates  FOXJ2,  HUWE 1 ,  NCOA5,  KRT 1 8,  TRA2B,  HNRNPM,  HNRNPA2B 1 , 
ZNF326  in  cell  culture. 

•  Isolation  of  an  in  vivo  NEK6  protein  complex  that  contains  two  kinases,  CSNKIAI  and  YESl,  which  are  also  in 
vitro  substrates  of  NEK6. 

•  Determination  that  NEK6  overexpression  in  tumors  leads  to  maintenance  of  immune-related  signaling  normally 
lost  with  castration,  as  well  as  squamous  transdifferentiation 

•  Demonstration  that  NEK6  substrates  FOXJ2,  CSNKIAI  and  YESl  mediate  a  subset  of  the  gene  expression 
changes  induced  by  NEK6  overexpression,  with  CSNKIAI  most  closely  related  to  the  squamous  differentiation 
phenotype 

•  Determination  that  NEK6  is  aberrantly  expressed  in  subset  of  human  prostate  cancer 

•  Optimization  of  a  genome-wide  barcoded  ORF  screen  of  genes  conferring  androgen  independent  proliferation  of 
LNCaP  cells  in  vitro 
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CONCLUSIONS 

We  have  performed  an  unbiased  in  vivo  functional  genomic  screen  to  identify  kinases  that  can  confer  androgen 
independence  in  a  model  of  androgen-dependent  prostate  tumor  formation,  and  discovered  NEK6  as  a  novel  mediator  of 
castration  resistance.  The  NEK6  gene  on  9q33.3  is  located  on  a  region  of  recurrent  copy  number  gain  in  prostate  cancer, 
and  it  is  overexpressed  in  a  subset  of  primary  human  prostate  cancers  and  most  patient-derived  prostate  cancer  cell  lines 
tested.  We  are  unable  to  correlate  protein  expression  to  amplification  of  the  locus  or  mRNA  overexpression  in  these 
samples,  so  it  is  unknown  whether  aberrant  NEK6  protein  expression  is  driven  by  genetic/transcriptional  changes  vs.  other 
translational  or  post-translational  determinants  of  protein  levels.  NEK6  plays  a  mechanistic  role  in  the  development  of 
castration  resistance  in  our  model,  and  turning  off  its  expression  in  xenograft  tumors  where  it  is  overexpressed  restores 
sensitivity  to  castration,  suggesting  that  its  continued  activity  is  required  for  tumor  maintenance  in  this  context.  The 
finding  that  NEK6  is  not  essential  for  mammalian  cells,  along  with  the  fact  that  mice  homozygous  for  a  targeted  null 
NEK6  allele  demonstrate  no  apparent  phenotype 

(https://www.infrafrontier.eu/sites/infrafrontier.eu/files/upload/public/deltagen/DELTAGEN  T518/)  except  for  promotion 
of  cardiac  hypertrophy  induced  by  transthoracic  aorta  constriction  (Bian  et  aL,  2014),  suggests  that  NEK6  could  be  safely 
therapeutically  targeted  in  prostate  cancer  and  other  cancer  types  where  it  is  a  relevant  driver  without  a  high  degree  of 
systemic  toxicity. 

The  role  of  NEK6  in  oncogenic  signaling  have  been  proposed  to  be  mediated  by  cell  cycle  progression/cyclin  B 
modulation  (Zhang  et  aL,  2014),  inhibition  of  p53-mediated  senescence  (Jee  et  aL,  2010),  phosphorylation  of  STAT3 
(Jeon  et  aL,  2010),  and  blocking  nuclear  translocation  of  SMAD4  to  antagonize  TGF(3  signaling  (Zuo  et  aL,  2014). 
However,  there  is  little  evidence  from  our  unbiased  studies  suggesting  any  of  these  pathways  as  the  predominant 
mechanism  of  NEK6-mediated  castration  resistance  in  our  model.  Previous  proteomic  work  (Vaz  Meirelles  et  aL,  2010) 
identifies  multiple  interacting  partners  and  substrates  that  implicate  NEK6  in  a  variety  of  signaling  pathways  including 
cell  cycle,  cytoskeleton  organization,  DNA  repair.  Notch  and  NFkB  pathways.  In  our  studies,  the  gene  expression 
signature  mediated  by  NEK6  overexpression  in  the  context  of  castration  implicates  differentiation  and  immune  processes 
in  conferring  castration  resistance  to  tumors.  Analysis  of  gene  expression  changes  induced  by  phosphomimic  forms  of 
newly  identified  substrates  FOXJ2,  NCOA5,  CKla  and  YESl  reveals  elements  of  NEK6  signaling  likely  mediated  by 
these  proteins.  In  particular,  the  gene  expression  profile  generated  through  expression  of  a  phosphomimic  form  of  CKla 
highly  correlates  with  the  NEK6  signature,  particularly  with  regards  to  markers  of  squamous  differentiation.  CKla  has 
previously  been  implicated  in  Wnt-|3-catenin  signaling  in  cancer  (Cruciat,  2014),  but  this  does  not  appear  to  be  the  operant 
mechanism  of  action  in  our  system  given  that  there  is  no  evidence  for  nuclear  (3-catenin  in  the  xenograft  tumors  by  IHC 
and  no  evidence  for  upregulation  of  a  P-catenin  signature  by  RNASeq/GSEA.  Recent  studies  suggest  CKla  as  a  potential 
therapeutic  target  in  acute  myeloid  leukemia  (Jaras  et  ai,  2014),  5q-  myelodysplastic  syndrome  (Schneider  et  al.,  2014), 
and  multiple  myeloma  (Hu  et  ai,  2015),  and  our  current  studies  would  suggest  NEK6/CKla  signaling  as  a  possible 
therapeutic  target  in  certain  forms  of  castration-resistant  prostate  cancer.  However,  given  that  the  substrates  identified 
here  are  not  sufficient  to  confer  castration  resistance,  there  are  likely  other  substrates  important  to  the  phenotype  that  we 
have  yet  to  characterize  (including,  for  example,  the  other  6  putative  substrates  identified  in  the  in  vitro  phosphoproteomic 
screen). 


Parallel  growth  factor  signaling  has  been  implicated  in  androgen  independence  in  several  model  systems,  but  this 
phenotype  depends  on  both  the  growth  factor  milieu  of  the  microenvironment  and  the  underlying  genetic  context  of  the 
cancer  cells  with  regards  to  their  behavior  in  this  milieu.  The  use  of  engineered  cell  lines  allows  control  over  the  genetic 
context,  and  the  use  of  xenografts  allows  testing  different  in  vivo  environments  where  the  limiting  nutrients  and  growth 
factors  are  unknown  and  thus  could  not  be  replicated  in  a  culture  environment.  While  NEK6  is  a  potent  mediator  of 
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castration  resistance  in  this  model,  it  is  reasonable  to  consider  that  other  genes  might  be  identified  to  confer  androgen 
independence  in  different  genetic  contexts  in  different  microenvironments. 

While  alterations  in  the  AR  gene  and  in  AR-mediated  signaling  have  already  been  demonstrated  to  play  a  role  in 
the  development  of  castration  resistance  (Sharifi,  2013),  it  is  apparent  that  resistance  to  hormonal  therapies  can  arise 
through  progression  to  states  where  the  androgen  receptor  is  no  longer  essential  for  survival,  collectively  termed  androgen 
pathway  independent  prostate  cancer  (APIPC)  (Nelson,  2012).  A  subset  of  these  cancers  express  markers  of 
neuroendocrine  differentiation,  and  transdifferentiation  to  a  squamous  phenotype  has  also  been  described  in  resistance  to 
hormonal  therapy  (Humphrey,  2012);  however,  the  molecular  mechanisms  underlying  transition  to  APIPC  are  not  fully 
characterized.  Specifically,  we  hypothesize  that  APIPC  represents  a  diversity  of  phenotypes  and  dependencies;  a 
comprehensive  understanding  of  these  mechanisms  would  identify  potential  therapeutic  targets  for  intervention,  and 
would  also  allow  for  biomarker  discovery  to  stratify  patients  for  likelihood  of  response  to  targeted  therapies  (Choudhury, 
Eeles,  et  al.,  2012).  Thus,  the  signaling  pathways  described  here  warrant  further  study  in  this  patient  population. 
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SUPPORTING  DATA 


Figure  1.  NEK6  confers  androgen-independent  tumor  formation  in  a  xenograft  model  of  androgen-dependent 
prostate  cancer.  A.  Tumor  formation  at  60  days  for  parental  LHSR-AR  cells  and  cells  expressing  NEK6  in 
female  and  castrated  mice.  B.  Tumor  formation  at  60  days  for  parental  LHMK-AR  cells  and  cells  expressing 
NEK6  in  female  mice.  C.  Inducible  expression  of  NEK6  in  vitro  at  48  hours  after  addition  of  doxycycline,  and 
waterfall  plot  of  change  in  tumor  volume  of  parental  LHSR-AR  cells  and  cells  with  inducible  NEK6  expression 
formed  in  male  mice,  30  days  after  castration  and  removal  of  testosterone  pellet. 

Figure  2.  A.  NEK6-mediated  androgen-independent  tumors  are  primarily  squamous  in  histology  and  AR 
negative.  Sections  of  tumors  derived  from  parental  LHSR-AR  cells  expressing  GFP  in  male  mice,  and  cells 
expressing  NEK6  in  female  and  castrated  mice  were  stained  with  AR  antibody  (brown).  B.  Sections  of  tumors 
derived  from  LHSR-AR  cells  overexpressing  NEK6  in  male  mice  with  implanted  testosterone  pellet  harvested 
prior  to  castration  (Day  0)  and  at  4,  8,  and  12  days  after  castration. 

Figure  3.  A.  Immunohistochemical  staining  at  20x  magnification  for  NEK6  in  prostate  cancer  tissue 
microarrays.  Low  (left  and  middle)  and  high  grade  (right)  cases  are  represented,  each  core  showing  tumor 
infiltrating  in  between  benign  glands  with  higher  expression  seen  in  the  tumor.  B.  Flow  diagram  for  assessment 
of  clinical  correlates  of  patients  with  high  vs.  low  NEK6  expression  in  tumors. 

Figure  4.  NEK6  is  overexpressed  in  several  prostate  cancer  cell  lines  compared  to  immortalized  (RWPE,  LH) 
and  transformed  (LHSR-AR)  prostate  epithelial  cells.  A.  Expression  of  NEK6  and  AR  in  prostate  cell  lines 
with  Hsp90  as  loading  control  B.  Expression  of  NEK6  in  CL-1,  PC-3,  and  DU  145  cells  with  doxycycline 
inducible  expression  of  2  shRNAs  targeting  NEK6  or  a  control  targeting  lacZ  in  the  presence  or  absence  of 
doxycycline  (CL-1  and  PC-3)  or  with  doxycycline  in  the  presence  or  absence  of  growth  factor  stimulation 
(DU145).  C.  Proliferation  curves  of  cells  cultured  in  the  presence  of  doxycycline  with  cells  in  60  mm  plates, 
split  and  replated  every  2-3  days  as  indicated.  Average  of  3  replicates  for  each  condition  is  shown. 

Figure  5.  A.  NEK6  can  phosphorylate  NCOA5  and  FOXJ2  in  vitro  at  the  sites  discovered  in  the 
phosphoproteomic  screen.  293T  cells  were  transfected  with  expression  constructs  for  wild-type  and  mutant  (S- 
to-D)  versions  of  NCOA5  or  FOXJ2  with  a  C-terminal  V5  tag  and  immunoprecipitated  with  anti-V5  antibody. 
Eluates  from  1/5  of  the  beads  were  assayed  by  V5  immunoblot;  the  remaining  4/5  was  subjected  to  on-bead  in 
vitro  kinase  assay  with  recombinant  active  GST-NEK6  (Sigma).  B+C.  Kinase  assay  performed  as  in  (A)  but 
with  wild-type  or  mutant  forms  of  FOXJ2  or  NCOA5  as  indicated:  +  indicates  a  mutation  is  present  at  that 
residue,  -  indicates  that  the  residue  is  wild-type  (lane  1  of  each  blot  represents  the  wild-type  protein).  5%  of  the 
input  assayed  by  V5  immunoblot  is  shown  in  the  left  panels;  the  kinase  assay  is  shown  in  the  right  panels. 
FOXJ2  has  10  S/T  residues  in  an  “acceptable”  NEK6  phosphorylation  motif;  the  six  previously  reported  in 
phosphosite.org  for  FOXJ2  or  homologous  to  a  previously  detected  phosphorylation  site  in  FOXJ3  were  tested 
in  this  assay.  NCOA5  has  5  S/T  residues  in  an  “acceptable”  NEK6  phosphorylation  motif;  however,  mutating 
Ser201  to  aspartic  acid  dramatically  decreased  exogenous  NCOA5  expression,  so  the  remaining  four  sites  were 
tested  in  this  assay. 

Figure  6.  A.  293T  cells  were  transfected  with  expression  constructs  for  the  indicated  kinases  with  a  C-terminal 
V5  tag  and  immunoprecipitated  with  anti-V5  antibody.  Eluates  from  1/4  of  the  beads  were  assayed  by  V5 
immunoblot;  the  remaining  3/4  was  subjected  to  on-bead  in  vitro  kinase  assay  with  recombinant  active  GST- 
NEK6  (Sigma).  B+C.  Kinase  assay  performed  as  in  (A)  but  with  wild-type  or  mutant  forms  of  CSNKl  A1  or 
YESl  as  indicated:  +  indicates  a  mutation  is  present  at  that  residue,  -  indicates  that  the  residue  is  wild-type  (lane 
1  of  each  blot  represents  the  wild-type  protein).  5%  of  the  input  assayed  by  V5  immunoblot  is  shown  in  the  left 
panels;  the  kinase  assay  is  shown  in  the  right  panels. 

Figure  7.  A.  Immunoblots  of  protein  lysates  of  xenograft  tumors  derived  from  LHSR-AR  cells  transduced  with 
doxycycline-inducible  constructs  for  the  expression  of  GFP  or  NEK6  in  male  mice,  with  doxycycline 
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maintained  in  the  diet  at  time  of  harvest  (+dox)  or  with  doxycycline  diet  removed  7  days  prior  to  harvest  (-dox). 
Tumors  were  harvested  from  non-castrated  mice  (day  0)  or  2  or  5  days  after  castration.  Two  exposures  of  the 
NEK6  immunoblot  are  shown.  B.  Venn  diagrams  (Venny  REF)  demonstrating  overlap  between  upregulated 
genes  (fold  change>1.5,  signal-to-noise>l)  and  downregulated  genes  (fold  change<0.67,  signal-to-noise<-l) 
with  NEK6  wild-type  vs.  kinase-dead  expression  in  LHSR-AR  cells  6  hours  after  growth  factor  stimulation  in 
vitro  assayed  by  Affymetrix  microarray  (green)  and  with  NEK6  overexpressing  vs.  control  LHSR-AR  xenograft 
tumors  2  days  (blue)  and  5  days  (yellow)  after  castration  assayed  by  RNASeq.  Overlapping  genes  are 
indicated.  C.  Gene  Set  Enrichment  Analysis  performed  on  the  gene  expression  signature  mediated  by  NEK6 
overexpression  at  day  5  after  castration  (GSEA  pre-ranked  based  on  ratio  of  classes  to  tumors  without  NEK6 
overexpression).  The  top  two  curated  gene  sets  (C2)  from  the  molecular  signatures  database  (MSigDB) 
correlated  with  the  NEK6  signature  are  shown. 

Figure  8.  A.  Gene  Set  Enrichment  Analysis  performed  on  the  gene  expression  signature  of  control  tumors  at  day 
5  after  castration  (GSEA  pre-ranked  based  on  ratio  of  classes  to  tumors  at  day  2  after  castration).  The  top  gene 
set  from  MSigDB  C2  correlated  with  genes  decreased  in  expression  after  castration  is  genes  involved  in 
interferon  alpha/beta  signaling  (http://www.reactome.org/cgi-bin/eventbrowser_st_id?ST_ID=REACT_25162); 
in  addition  two  prostate  cancer  related  gene  sets  are  also  highly  correlated  with  this  signature  (ranked  #’s  9  and 
12):  genes  up-regulated  in  prostate  tumors  developed  by  transgenic  mice  overexpressing  VAV3  (Liu  et  ah, 
2008)  and  genes  up-regulated  in  prostate  cancer  samples  from  African-American  patients  compared  to  those 
from  the  European-American  patients  (Wallace  et  al.,  2008).  B.  Genes  downregulated  in  the  control  (-dox) 
tumors  after  castration  (fold  change  <-1.5,  signal-to-noise  <-l  in  the  comparison  of  tumors  harvested  at  day  5 
and  day  2)  were  plotted  against  the  NEK6  signature  at  day  5  using  GSEA  pre-ranked  (left)  and  compared  with 
genes  upregulated  by  NEK6  at  day  5  (fold  change  >1.5,  signal-to-noise  >1  in  the  comparison  of  +dox  to  -dox 
tumors)  by  Venn  diagram  (right).  The  GO  terms  most  highly  enriched  in  this  overlap  are  cytokine-mediated  and 
type  I  interferon  signaling  pathways  (p=2xl0‘^).  C.  RNASeq  was  performed  from  3  tumors  inducibly 
expressing  a  phosphomimic  form  of  CKla  (S96D-I-T129D-I-T287D)  with  continued  doxycycline  expression,  and 
3  matched  tumors  s/p  doxycycline  withdrawal  5  days  after  castration.  Left:  The  correlation  of  genes  upregulated 
by  CKla(Dx3)  with  the  NEK6  signature  at  day  5  by  GSEA  pre-ranked;  middle:  overlap  of  the  upregulated 
genes;  right:  the  correlation  of  the  CKla(Dx3)  signature  with  the  Bosco  Epithelial  Differentiation  Module 
(right). 


Supplemental  Figure  1.  A.  Immunoblots  using  Epitomics  3789-1  (left)  and  Santa  Cruz  H-50  (right)  anti-NEK6 
primary  antibodies  on  protein  lysates  of  LHSR-AR  cells  transduced  with  a  doxycycline-inducible  construct  for 
the  expression  of  NEK6  with  (-I-)  or  without  (-)  the  addition  of  doxycycline  in  the  media  for  48  hours.  B. 
Immunohistochemistry  of  xenograft  tumors  derived  from  LHSR-AR  cells  transduced  with  doxycycline- 
inducible  constructs  for  the  expression  of  NEK6  in  male  mice,  with  doxycycline  maintained  in  the  diet  at  time 
of  harvest  (+doxycycline)  or  with  doxycycline  diet  removed  7  days  prior  to  harvest  (-doxycycline).  Tumors 
were  harvested  5  days  after  castration.  C.  Percentage  of  NEK6  positive  cells  (defined  as  having  cytoplasmic 
NEK6  staining  intensity  greater  than  two  standard  deviations  above  the  median)  plotted  against  mean 
cytoplasmic  NEK6  staining  intensity  (normalized  as  a  percentile  rank  in  each  TMA)  of  cells  within 
histologically  benign-appearing  glandular  tissue  (left)  and  tumor  tissue  (right)  from  radical  prostatectomy 
specimens  of  patients  with  prostate  cancer.  Each  point  represents  composite  data  from  all  benign  or  tumor 
cores  from  an  individual  patient.  The  threshold  for  positivity  was  set  as  the  discontinuity  in  the  distribution  of 
%  positive  cells  (at  12%)  as  indicated  by  the  dashed  red  lines. 

Supplemental  Figure  2.  A.  Luciferase  activity  detected  in  LNCaP  cells  transiently  transfected  with  an  AR 
reporter  either  alone  (top)  or  in  combination  with  an  expression  plasmid  for  NEK6  (bottom)  incubated  with 
concentrations  of  the  synthetic  androgen  R1 881  indicated.  B.  Doxycycline-inducible  expression  of  NEK6 
wild- type,  kinase  dead  (K74M/K75M),  AKTl  and  RAFl  both  untagged  and  with  C-terminal  V5  tag.  Untagged 
versions  were  used  for  experiments  described.  C.  Inducible  expression  of  NEK6  does  not  increase  expression 
of  AR  targets  PSA  or  TMPRSS2  in  LHSR-AR  cells.  Expression  of  TMPRSS2  and  PSA  were  measured  by 
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qPCR  in  the  absence  and  presence  of  doxycycline  to  induce  transgene  expression  and  in  the  presence  and 
absence  of  R1881  as  indicated,  with  expression  normalized  to  cells  transduced  with  NEK6  K74M/K75M  in  the 
absence  of  doxycycline  and  R1881.  D.  NEK6  expression  has  neither  a  positive  or  negative  effect  on  AR 
signaling  as  measured  through  published  AR  signatures.  Gene  expression  changes  conferred  by  inducible 
expression  of  wild-type  NEK6  vs.  kinase  dead  NEK6  six  hours  after  growth  factor  stimulation  were  assayed  in 
3  biological  replicates,  and  GSEA  was  used  to  assess  enrichment  of  signatures  positively  correlated  with  AR 
activity  in  two  data  sets  from  the  literature  (Hieronymus  et  al.,  2006;  Mendiratta  et  ah,  2009,  Sharma  et  ah, 
2013). 

Supplemental  Figure  3.  A.  Immunohistochemistry  of  xenograft  tumors  derived  from  LHSR-AR  cells 
transduced  with  doxycycline- inducible  constructs  for  the  expression  of  NEK6  in  male  mice,  with  doxycycline 
maintained  in  the  diet  at  time  of  harvest  (+doxycycline)  or  with  doxycycline  diet  removed  7  days  prior  to 
harvest  (-doxycycline).  Tumors  were  harvested  from  non-castrated  mice  (day  0)  or  2  or  5  days  after  castration. 
B.  LHSR-AR  cells  expressing  GFP  or  wild-type  or  mutant  forms  of  NEK6  under  a  doxycycline  inducible 
promoter,  in  the  presence  of  doxycycline  and  absence  of  growth  factor  stimulation  (lanes  1-4),  with  growth 
factor  stimulation  in  the  absence  of  doxycycline  (lanes  5-9)  and  with  growth  factor  stimulation  in  the  presence 
of  doxycycline  (lanes  10-14).  Immunoblots  for  phosphorylated  forms  of  p70S6K  at  Thr412  and  STAT3  at 
Ser727  demonstrate  increases  with  growth  factor  stimulation  (compare  lanes  1-4  with  lanes  10-14)  but  no 
increase  in  the  presence  of  wild- type  NEK6  (lanes  2,  12)  as  compared  to  kinase-dead  NEK6  (lanes  1,  11). 

Supplemental  Figure  4.  NEK6  overexpression  does  not  lead  to  promotion  of  cell  cycle  progression  or 
antagonism  of  the  p53  pathway  in  LHSR-AR  cells.  A.  Cell  cycle  profiles  of  cells  with  and  without  NEK6 
overexpression.  LHSR-AR  cells  transduced  with  doxycycline-inducible  NEK6  were  cultured  in  the  presence 
(bottom  profiles)  or  absence  (top  profiles)  of  doxycycline,  starved  from  growth  factors  for  24  hours,  and 
released  into  growth  factor-containing  media  for  the  times  indicated,  then  harvested  and  fixed  for  propidium 
iodide  staining  in  comparison  to  asynchronously  cycling  cells  B.  Proliferation  curves  of  LHSR-AR 
constitutively  expressing  lacZ  or  NEK6,  average  cell  counts  from  3  plates  collected  at  the  indicated  time  points 
plotted  compared  to  previous  time  point  with  standard  deviations.  C.  p53  pathway  is  inactive  in  LHSR-AR 
cells,  and  NEK6  expression  does  not  rescue  from  cell  death  mediated  by  etoposide.  LHSR-AR  cells  with 
constitutive  expression  (pLX304-)  of  lacZ  vs.  NEK6  or  doxycycline-inducible  expression  (pTRIPz-)  of  NEK6 
kinase  dead  (kd)  vs.  wild-type  (wt)  were  exposed  to  etoposide  at  concentrations  indicated  vs.  DMSO  as  vehicle 
control.  Attached  and  floating  cells  were  harvested  and  combined  for  immunoblotting. 

Supplemental  Figure  5.  RNASeq  was  performed  from  3  tumors  inducibly  expressing  phosphomimic  forms  of 
A.  FOXJ2  (S8D+T23D+S254D)  B.  YESl  (T295D+S524D)  and  C.  NCOA5  (S21D+S96D+S151D)  with 
continued  doxycycline  expression,  and  3  matched  tumors  for  each  s/p  doxycycline  withdrawal  5  days  after 
castration.  The  genes  upregulated  by  FOXJ2(Dx3)  and  YESl(Dx2),  but  not  NCOA5(Dx3),  were  highly 
correlated  with  the  NEK6  signature  at  day  5  by  GSEA  pre-ranked  (left),  with  overlapping  upregulated  genes 
with  NEK6  shown  in  the  Venn  diagrams  (right).  Blue  represents  differentiation  process  genes,  red  represents 
immune  process  genes,  purple  represents  genes  in  both  categories,  and  black  represents  genes  in  neither 
category. 

Supplemental  Figure  6.  A-i-B.  293T  cells  were  transfected  with  expression  constructs  for  the  indicated  proteins 
with  a  C-terminal  V5  tag  and  immunoprecipitated  with  anti-V5  antibody.  Eluates  from  1/4  of  the  beads  were 
assayed  by  V5  immunoblot;  the  remaining  3/4  was  subjected  to  on-bead  in  vitro  kinase  assay  with  recombinant 
active  GST-NEK6  (Sigma). 
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Supplemental  Table  1. 


AR  Nuclear  % 
score  (0,1, 2,3,4)" 

AR  Nuclear 
intensity^ 

AR% 

Cytoplasm  score 
(0,1, 2, 3,4)" 

AR  Cytoplasm 
intensity^ 

notes 

Male  mice 

(positive 

control) 

GFP  1  (one 
piece) 

4 

m  s 

3 

w  m 

9.11.12 

GFP  1  (second 
piece) 

3 

w  m  s 

1 

w 

GFP  2 

4 

m  s 

1 

w 

GFP3 

4 

s 

4 

m  s 

Female 

mice 

NEK6  1  (one 
piece) 

2 

m  s 

2 

w  m 

NEK6  2  (few 
pieces) 

1 

w 

1 

w 

NEK6  3 

trace 

most  nuclei  are 
neg;  occasional 
nuclei  show  weak 
staining  (<1-2%) 

trace 

most  cytoplasm 
are  neg; 
occasional  cells 
show  weak 
staining  (<1-2%) 

surface  mouse 
skin  is  neg  for  AR 

AKT1 

2 

m  s 

3 

w  m  s 

CCL2 

0  (very  small 
sample) 

0 

0  (very  small 
sample) 

0 

scant  epithelial 
cells 

ERBB2 

2 

w  m  s 

2 

w  m 

KRASV12 1 
(one  piece) 

0 

0 

KRASV12  2(2 
other  pieces) 

1 

m  s 

2 

w  m 

MEKDD 

1 

w 

4 

w  m  s 

PIM1 

<5% 

w  (rare  strong,  3- 
4  cells) 

2 

m  s 

RAF1 

1 

w  m 

4 

w  m 

two  pieces  show 
this  staining 
pattern:  third 
piece  has  weak 
cyto  stain  in 
<25%,  and 
nuclear  stain 
weak  in  <25%. 

Castrated 

mice 

NEK6  castr 

most  are  neg; 
however,  one 
tissue  fragment 
shows  up  to  3+ 
Nuc  score  in 
-25%  of  nuclei 

0-3 

^AR  %  score 

0,  1(1 -25%), 2(25-50%),  3(50-75%),  4(75-100%) 

*  Intensity:  \n- 

=Weak,  m=Moderate,  s=Strong  (bold  indicates  predominant  intensity) 
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Supplemental  Table  2. 


Gene  Name 

Site 

q  enriched  (wt 
vs.  kd) 

q  enriched  (wt 
induced  vs. 
uninduced) 

LIMCHI 

S303 

9.73E-64 

1.72E-102 

DROSHA 

S357 

7.81E-32 

2.87E-68 

FOXJ2 

S8 

4.12E-18 

2.22E-06 

F0X03 

S7 

1.06E-17 

4.06E-54 

HUWEl 

S2595 

7.84E-13 

1.09E-09 

EPPKl 

S1529 

4.20E-12 

2.03E-06 

COIL 

S487 

6.79E-12 

1.13E-13 

SRGAPl 

S932 

6.79E-12 

1.21E-25 

SLC2A12 

S244 

6.95E-11 

2.89E-05 

TRPSl 

S843 

2.12E-09 

3.99E-13 

MTXl 

S9 

4.78E-06 

3.80E-03 

OGFR 

S349 

9.83E-06 

9.02E-07 

INTS3 

S502 

5.00E-05 

2.61E-06 

TRA2B 

S239 

7.67E-05 

3.27E-03 

LM07 

S926 

9.02E-05 

5.17E-02 

SATB2 

S20 

9.20E-05 

5.42E-02 

ATM 

T1885 

1.32E-04 

1.83E-02 

PLEKHA6 

S313 

3.42E-04 

2.86E-06 

LM07 

S895 

3.62E-04 

2.13E-01 

PAK6 

S246 

4.80E-04 

8.11E-07 

HNRNPM 

S633 

6.05E-04 

3.70E-04 

MLLT3 

S302 

2.20E-03 

9.42E-02 

HNRNPA2B1 

S324 

2.41E-03 

5.43E-03 

NCOA5 

S96 

2.52E-03 

1.18E-01 

SETX 

T2474 

2.53E-03 

1.92E-01 

BCL6 

S466 

2.65E-03 

2.40E-01 

ATXNl 

S811 

2.94E-03 

1.27E-03 

CDKN2AIP 

S151 

3.23E-03 

1.47E-01 

FAM21C 

S288 

3.49E-03 

1.82E-01 

EPS8L1 

T305 

5.41E-03 

3.80E-02 

FOXAl 

S307 

5.42E-03 

l.OOE-01 

KLF4 

T316 

6.40E-03 

6.99E-05 

LM07 

S1593 

6.56E-03 

1.32E-02 

ZNF326 

S131 

8.20E-03 

1.57E-01 

PLEKHG6 

S645 

2.05E-02 

1.29E-02 

RIPK3 

S316 

3.04E-02 

1.70E-03 

ZNF326 

S106 

3.98E-02 

9.73E-09 

ERCC5 

S156 

4.94E-02 

1.03E-01 

ERCC5 

S157 

4.94E-02 

1.03E-01 

LIGl 

S66 

5.02E-02 

9.97E-02 

MYOF 

S193 

5.29E-02 

9.46E-04 

ZDHHC18 

S19 

6.19E-02 

1.13E-07 

PBRMl 

S353 

6.53E-02 

1.45E-09 

modifiedsequencejocalization 

_SWSTATS(ph)PLGGERPFR_7 

_NTDSWAPPLEIVNHRS(ph)PS(ph)REK_18 

_{ac)ASDLESS{ph)LTSIDWLPQLTLR_7 

_{ac)AEAPAS(ph)PAPLSPLEVELDPEFEPQSRPR_6 

_LLGPSAAADILQLSSS{ph)LPLQSR_16 

_QVS(ph)ARDLFR_3 

_KIDS{ph)PPIRR_4 

_LLELTS{ph)SYSPDVSDYKEGR_7 

_LRALS{ph)DTTEELTVIK_5 

_TLRDS(ph)PNVEAAHLARPIYGLAVETK_5 

_{ac)MLLGGPPRS{ph)PR_9 

_S{ph)VEPQDAGPLER_l 

_FPEFCSSPS{ph)PPVEVK_9 

_S(ph)YRGGGGGGGGWR_l 

_GISS(ph)LPR_4 

_SGS{ph)PDVKGPPPVK_3 

_STT(ph)PANLDSESEHFFR_3 

_KSS{ph)MNQLQQWVNLRR_3 

_VSAS(ph)LPR_4 

_HGSEEARPQSCLVGSATGRPGGEGS(ph)PS(ph)PK_25 

_GNFGGS(ph)FAGSFGGAGGHAPGVAR_6 

_KKS(ph)SSEALFK_3 

_SGNFGGS{ph)RNMGGPYGGGNYGPGGSGGSGGYGGR_7 

_DLRDS(ph)RDFR_5 

_SLT{ph)HPPTIAPEGSRPQGGLPSSKLDSGFAK_3 

_SSSESHS{ph)PLYMHPPK_7 

JCIEGRS(ph)NVGK_7 

_VIEGKNS(ph)SAVEQDHAK_7 

_S{ph)RPTS(ph)FADELAAR_5 

_AAGEGLLT{ph)LR_8 

_KDPSGASNPSADS(ph)PLHR_13 

_TT(ph)PTLGLEEVLSSR_2 

_SHS(ph)PSASQSGSQLR_3 

_NQGGSS(ph)WEAPYSR_6 

_S{ph)APELPEGILK_l 

_RFS{ph)IPESGQGGTEMDGFRR_3 

_FGGS(ph)YGGRFESSYR_4 

_ENDLYVLPPLQEEEKHS(ph)S(ph)EEEDEKEWQER_17 

_ENDLYVLPPLQEEEKHS(ph)S(ph)EEEDEKEWQER_18 

_VLGS{ph)EGEEEDEALS{ph)PAK_4 

_RMLS(ph)NKPQDFQIR_4 

_{ac)MKDCEYQQISPGAAPLPAS(ph)PGAR_19 

_LSAITM(ox)ALQYGS{ph)ES(ph)EEDAALAAAR_12 
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EXPH5 

S1444 

7.11E-02 

3.82E-02 

SIPA1L3 

S158 

7.11E-02 

5.33E-07 

ATXNl 

S238 

7.85E-02 

9.85E-02 

KRT18 

S323 

8.16E-02 

1.44E-01 

IRF2BP1 

S453 

9.61E-02 

4.17E-02 

DLG3 

Y673 

l.OlE-01 

1.62E-01 

ATXNl 

S775 

1.02E-01 

7.46E-03 

EXPH5 

S341 

1.12E-01 

3.01E-02 

RFX2 

S28 

1.27E-01 

2.98E-02 

PCYTIB 

S315 

1.59E-01 

3.42E-02 

EPS8L1 

T202 

1.61E-01 

2.35E-01 

HIVEP2 

S2300 

1.81E-01 

3.82E-02 

ZNF608 

S964 

1.96E-01 

4.92E-05 

KLF3 

S71 

2.08E-01 

2.39E-23 

ACLY 

S481 

2.24E-01 

2.45E-04 

ZNF608 

S1453 

2.31E-01 

4.83E-07 

RSS(ph)WECTGSGR_3 

SKDVEFQDGWPRS{ph)PGR_13 

APGLITPGS(ph)PPPAQQNQYVHIS{ph)SSPQNTGR_9 

NLKASLENS{ph)LREVEAR_9 

NVAEALGHSPKDPGGGGGPVRAGGAS{ph)PAASSTAQPPTQHR_26 

RDNEVDGQDY(ph)HFVVSR_10 

WS(ph)APESR_2 

S{ph)LHFPATTQSK_l 

.{ac)MQNSEGGADSPASVALRPSAAAPPVPAS(ph)PQR_28 

M{ox)LQALS(ph)PK_6 

AVIST{ph)VER_5 

RGPHALQSSGPPSTPS(ph)SPR_17 

SKASS{ph)PSDIISSKDSVVK_5 

S{ph)SPPSAGNSPSSLKFPSSHRR_2 

KAKPAMPQDSVPS(ph)PR_13 

DRHS{ph)PFGQR_4 
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Supplemental  Table  3. 


GeneNames 

log2  samplexontrol  median 

totalintensity 

numSpectra 

Unique  peptides 

NEK6 

2.603 

2.48E+12 

118 

16 

ACTR2 

2.419 

2.44E+12 

42 

17 

ACTRIA 

2.426 

2.41E+12 

8 

4 

MYH9 

2.651 

1.49E+12 

291 

61 

MYL6 

2.724 

1.28E+12 

42 

8 

UBA52 

1.995 

2.46E+11 

8 

4 

SFN 

1.987 

1.79E+11 

22 

8 

FRASl 

4.813 

1.53E+11 

2 

2 

COROIA 

3.79 

9.88E+10 

3 

2 

GSN 

2.555 

7.40E+10 

41 

14 

EEFlAl 

2.025 

6.04E+10 

45 

13 

S100A8 

1.85 

5.65E+10 

27 

12 

TUBB 

2.159 

4.58E+10 

14 

3 

S100A9 

1.325 

4.51E+10 

21 

10 

MYHIO 

2.501 

3.80E+10 

34 

12 

TUBB2C 

2.115 

3.18E+10 

20 

4 

MY06 

2.184 

3.05E+10 

41 

20 

ARPC2 

2.378 

2.97E+10 

40 

17 

MYH14 

1.85 

2.64E+10 

16 

7 

UTP20 

6.005 

2.09E+10 

4 

2 

ARPC4 

2.567 

2.05E+10 

22 

7 

DBT 

1.389 

2.03E+10 

21 

9 

FLNB 

1.724 

1.80E+10 

53 

29 

TMODl 

3.374 

1.74E+10 

8 

2 

FLNA 

2.156 

1.70E+10 

55 

22 

TMOD3 

2.201 

1.64E+10 

27 

11 

IQGAPl 

1.641 

1.51E+10 

27 

15 

EIF4A1 

1.931 

1.49E+10 

21 

10 

MYOIE 

2.371 

1.41E+10 

19 

8 

KIF5B 

1.293 

1.33E+10 

7 

6 

HSPA8 

1.543 

1.32E+10 

41 

13 

CNGBl 

3.565 

1.26E+10 

3 

2 

CAPZB 

1.846 

1.07E+10 

21 

13 

ACTN4 

2.523 

1.02E+10 

14 

7 

MYL9 

1.78 

8.66E+09 

12 

4 

PRDX4 

2.438 

7.98E+09 

9 

2 

CDC37 

2.483 

7.74E+09 

24 

13 

S100A4 

3.195 

7.18E+09 

8 

4 

TUBB2A 

1.997 

7.12E+09 

8 

2 

MY018A 

2.216 

6.92E+09 

19 

17 

SLC25A3 

3.05 

6.71E+09 

10 

5 

S100A2 

3.209 

6.39E+09 

11 

6 

MYOID 

2.363 

6.20E+09 

19 

7 

TPIl 

1.311 

5.70E+09 

15 

8 

PSMDl 

1.932 

5.45E+09 

3 

3 

ANXA2 

1.802 

5.22E+09 

14 

9 

MYOIB 

2.951 

5.07E+09 

10 

5 

ACTCl 

2.188 

5.03E+09 

12 

2 

CALMl 

2.67 

4.93E+09 

7 

4 

HSP90AA1 

2.25 

4.81E+09 

14 

10 

TPMl 

3.045 

4.78E+09 

6 

5 

NPMl 

2.247 

4.67E+09 

6 

5 

ARPC5 

1.533 

4.55E+09 

12 

6 

PRKDC 

1.331 

4.53E+09 

16 

12 

SERPINB5 

1.815 

4.48E+09 

17 

11 

EEF2 

1.533 

4.29E+09 

20 

15 

TPM4 

1.951 

3.89E+09 

7 

5 

PGAM5 

1.324 

3.84E+09 

3 

2 

YWHAZ 

2.909 

3.55E+09 

6 

1 

CRIPl 

1.698 

3.41E+09 

5 

3 

RAN 

1.751 

3.30E+09 

11 

6 

THOC3 

2.879 

3.23E+09 

4 

2 

PLS3 

1.543 

3.19E+09 

16 

9 

NEK9 

2.06 

2.91E+09 

14 

9 

TPM4 

2.075 

2.87E+09 

5 

3 

TUFM 

1.592 

2.82E+09 

12 

6 

DSTN 

2.356 

2.68E+09 

3 

2 

RUVBL2 

1.737 

2.65E+09 

10 

8 

ANXAl 

1.329 

2.64E+09 

15 

7 

SPRRIB 

1.522 

2.49E+09 

9 

6 

YWHAG 

2.889 

2.48E+09 

4 

3 
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SPTBN2 

2.007 

2.48E+09 

7 

7 

ALDOC 

1.406 

2.39E+09 

3 

2 

DOCK? 

2.047 

2.32E+09 

8 

6 

DNAJA2 

1.743 

2.23E+09 

8 

5 

RUVBLl 

2.108 

2.17E+09 

14 

6 

GNBl 

2.124 

2.02E+09 

6 

3 

PCBP2 

1.573 

1.94E+09 

6 

3 

CRIP2 

2.268 

1.91E+09 

8 

3 

RAB7A 

1.509 

1.88E+09 

14 

6 

SERPINBl 

1.626 

1.85E+09 

6 

4 

TCPl 

2.157 

1.82E+09 

11 

8 

ARPCIA 

2.279 

1.79E+09 

6 

6 

LIMAl 

2.646 

1.77E+09 

7 

6 

GNB2 

2.532 

1.73E+09 

5 

3 

SLC2A1 

3.03 

1.72E+09 

2 

2 

CLTA 

2.27 

1.70E+09 

5 

4 

VCP 

1.472 

1.66E+09 

12 

6 

HNRNPK 

1.417 

1.57E+09 

5 

4 

SLC25A6 

1.831 

1.44E+09 

9 

3 

RABllB 

1.291 

1.40E+09 

8 

5 

SEC61A1 

2.093 

1.39E+09 

6 

4 

ZNF185 

2.191 

1.37E+09 

5 

5 

PPIA 

1.374 

1.33E+09 

17 

10 

NPEPPS 

1.421 

1.32E+09 

12 

10 

AIFIL 

2.351 

1.30E+09 

3 

3 

GNB2L1 

1.616 

1.29E+09 

10 

9 

SETX 

2.13 

1.24E+09 

4 

2 

PHGDH 

1.6 

1.24E+09 

6 

6 

KIF21B 

2.038 

1.21E+09 

4 

3 

AHNAK 

1.377 

1.21E+09 

7 

6 

DNAJA4 

2.71 

1.20E+09 

2 

1 

CCT3 

1.365 

1.20E+09 

9 

6 

MAPKl 

1.597 

1.19E+09 

3 

3 

PPPICA 

2.184 

1.15E+09 

4 

4 

AHCY 

1.393 

1.14E+09 

7 

4 

PDCD6IP 

1.369 

1.12E+09 

17 

11 

STT3A 

2.571 

l.llE+09 

6 

4 

TEX15 

2.183 

l.lOE+09 

3 

3 

SQRDL 

1.705 

1.08E+09 

4 

3 

CALML3 

2.365 

1.04E+09 

7 

4 

NCL 

1.393 

1.04E+09 

9 

5 

MY05A 

2.628 

1.02E+09 

5 

3 

SLC25A11 

1.729 

9.75E+08 

7 

3 

HP1BP3 

1.362 

9.71E+08 

7 

6 

TUBB6 

1.576 

9.46E+08 

9 

3 

YMEILI 

1.565 

9.32E+08 

2 

2 

CPNE3 

1.498 

8.91E+08 

5 

4 

ATL3 

1.423 

8.90E+08 

2 

2 

PHB 

1.577 

8.85E+08 

5 

3 

DDOST 

1.817 

8.72E+08 

4 

3 

CAMK2D 

1.354 

8.71E+08 

9 

6 

SERPINB8 

1.446 

8.54E+08 

3 

3 

HARS 

2.995 

8.41E+08 

3 

2 

COPB2 

2.197 

7.95E+08 

6 

5 

POTEKP 

1.725 

7.90E+08 

3 

1 

TMEDIO 

1.515 

7.88E+08 

3 

3 

SIOOAIO 

1.918 

7.68E+08 

2 

1 

KRAS 

1.285 

7.63E+08 

3 

2 

NARS 

1.771 

7.54E+08 

4 

3 

H2AFJ 

1.726 

7.44E+08 

2 

1 

MYL12B 

2.834 

7.31E+08 

2 

1 

CTTN 

1.64 

6.95E+08 

7 

7 

KPNBl 

1.733 

6.91E+08 

3 

3 

LADl 

1.318 

6.69E+08 

3 

2 

YWHAH 

2.803 

6.58E+08 

4 

3 

HSPHl 

2.016 

6.27E+08 

8 

6 

SERPINB13 

1.343 

5.93E+08 

3 

3 

SSR4 

1.925 

5.92E+08 

2 

1 

COPBl 

2.574 

5.86E+08 

5 

3 

MY05C 

2.939 

5.66E+08 

3 

3 

PDLIMl 

1.423 

5.53E+08 

4 

4 

PDLIM5 

1.733 

5.50E+08 

4 

2 

PRDX6 

2.199 

5.49E+08 

2 

2 

29 


MCM5 

1.941 

5.49E+08 

2 

2 

TPMl 

3.063 

5.31E+08 

4 

2 

XPOl 

1.629 

5.31E+08 

3 

3 

C22orf28 

1.391 

5.18E+08 

5 

5 

SDR16C5 

1.795 

5.00E+08 

6 

3 

PSMC6 

1.433 

4.99E+08 

4 

3 

NEK7 

1.287 

4.97E+08 

4 

1 

GSTPl 

1.323 

4.83E+08 

6 

3 

ATP50 

1.373 

4.81E+08 

6 

3 

HSPA9 

1.453 

4.76E+08 

6 

4 

CCT4 

1.418 

4.73E+08 

5 

2 

HIST1H2BB 

1.991 

4.71E+08 

2 

1 

CDK5RAP2 

1.542 

4.62E+08 

2 

2 

ALDHIAI 

1.93 

4.50E+08 

2 

2 

TUBAIC 

1.418 

4.41E+08 

3 

1 

DHCR7 

1.404 

4.24E+08 

3 

2 

EIF4A3 

1.442 

4.16E+08 

3 

2 

NAPILI 

1.771 

4.07E+08 

3 

1 

GMFB 

1.831 

4.05E+08 

4 

3 

CYFIP2 

1.359 

4.03E+08 

5 

5 

PPP2R2A 

1.373 

4.01E+08 

4 

3 

CDKALl 

2.115 

3.99E+08 

2 

2 

RAB5A 

1.3 

3.99E+08 

4 

2 

ATP5C1 

1.371 

3.90E+08 

4 

3 

CKB 

1.775 

3.89E+08 

4 

4 

PSMC5 

1.409 

3.75E+08 

3 

3 

CAPZAl 

1.302 

3.63E+08 

3 

2 

SLC25A13 

1.574 

3.60E+08 

3 

1 

NME2 

1.266 

3.60E+08 

4 

2 

H2AFZ 

1.392 

3.49E+08 

2 

2 

PSMCl 

2.106 

3.48E+08 

4 

3 

HECTDl 

1.35 

3.46E+08 

4 

4 

PPIB 

1.443 

3.44E+08 

4 

4 

CTNNDl 

1.785 

3.42E+08 

2 

2 

SAMD9 

2.206 

3.40E+08 

5 

5 

IKBKB 

1.889 

3.38E+08 

2 

2 

ADH5 

1.927 

3.30E+08 

2 

2 

UACA 

1.533 

3.29E+08 

4 

4 

ATPlAl 

1.685 

3.17E+08 

2 

1 

SlOOAll 

1.47 

3.13E+08 

5 

2 

C5orf38 

2.752 

3.10E+08 

2 

2 

ATAD3B 

1.253 

3.10E+08 

2 

2 

CEP350 

2.936 

3.09E+08 

2 

2 

RAB14 

1.999 

3.09E+08 

3 

3 

RPS6KA3 

2.241 

3.06E+08 

3 

3 

lARS 

2.182 

3.06E+08 

2 

2 

DCTNl 

1.341 

2.96E+08 

2 

2 

YESl 

2.071 

2.92E+08 

3 

3 

TAGLN2 

1.69 

2.90E+08 

3 

3 

MAP2K2 

1.646 

2.86E+08 

2 

2 

DDBl 

1.72 

2.83E+08 

3 

3 

XDH 

1.432 

2.78E+08 

3 

3 

PCFll 

1.932 

2.74E+08 

2 

2 

NAMPT 

1.75 

2.74E+08 

3 

3 

CSNKIAI 

2.1 

2.55E+08 

2 

1 

MT-C02 

1.385 

2.52E+08 

2 

2 

SVIL 

2.039 

2.51E+08 

2 

2 

ESYT2 

1.444 

2.51E+08 

3 

3 

PSMD14 

1.693 

2.50E+08 

2 

2 

SFRS3 

2.106 

2.42E+08 

2 

1 

YWHAE 

2.248 

2.38E+08 

2 

2 

SEPT2 

1.854 

2.37E+08 

3 

2 

RPS6KA1 

2.005 

2.25E+08 

3 

3 

COROIB 

1.378 

2.24E+08 

4 

2 

RABIB 

1.89 

2.17E+08 

3 

2 

SNRPD3 

1.277 

2.15E+08 

2 

2 

HSPA2 

1.325 

2.07E+08 

2 

2 

VPS35 

2.32 

1.96E+08 

2 

2 

VPS29 

1.264 

1.95E+08 

2 

2 

AHCYLl 

1.289 

1.93E+08 

2 

2 

DDX3X 

2.104 

1.91E+08 

2 

2 

ACTR3 

2.319 

1.86E+08 

2 

2 

TWFl 

1.427 

1.85E+08 

3 

2 
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FNl 

1.328 

1.85E+08 

3 

3 

UNC45A 
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GO  terms  NEK6  wt  vs.  NEK6  kd  +growtti  faetor  stimulation  in  vitro 


E 


GOID 

GOTERM 

Level 

RVALUE 

Value 

GeneHits 

G0:0051240 

positive  regulation  of  multicellular  organismal  process 

3 

0.008297 

3.42 

6 

G0:0008015 

blood  circulation 

5 

0.009059 

3.92 

5 

G0:0003013 

circulatory  system  process 

4 

0.009197 

3.91 

5 

G0:0065008 

regulation  of  biological  quality 

3 

0.010748 

1.73 

19 

G0:0048878 

chemical  homeostasis 

5 

0.012744 

2.56 

8 

G0:0009968 

negative  regulation  of  signal  transduction 

4 

0.016305 

2.95 

6 

G0:0023057 

negative  regulation  of  signaling 

3 

0.022366 

2.74 

6 

G0:0010648 

negative  regulation  of  cell  communication 

4 

0.022577 

2.74 

6 

G0:0010817 

regulation  of  hormone  levels 

4 

0.023326 

3.08 

5 

G0:0051239 

regulation  of  multicellular  organismal  process 

3 

0.023671 

1.89 

12 

GO  terms  interseetion  of  Castration  Down  (Day  2  to  5)  and  NEK6  s/p  eastration  Day  5  Up 

F 

GOID 

GOTERM 

Level 

RVALUE 

Value 

GeneHits 

G0:0045087 

innate  immune  response 

4 

0.000041 

7.31 

7 

G0:0052548 

regulation  of  endopeptidase  activity 

7 

0.000092 

10.84 

5 

G0:0052547 

regulation  of  peptidase  activity 

6 

0.000108 

10.48 

5 

G0:0051336 

regulation  of  hydrolase  activity 

5 

0.000419 

5.01 

7 

G0:0006955 

immune  response 

3 

0.000871 

3.87 

8 

G0:0006952 

defense  response 

4 

0.000988 

3.79 

8 

G0:0009607 

response  to  biotic  stimulus 

3 

0.001305 

4.87 

6 

G0:0006950 

response  to  stress 

3 

0.002619 

2.29 

13 

G0:0042221 

response  to  chemical  stimulus 

3 

0.003219 

2.35 

12 

G0:0051707 

response  to  other  organism 

3 

0.003466 

4.85 

5 

GO  terms  interseetion  of  Castration  Up  (Day  2  to  5)  and  NEK6  s/p  eastration  Day  5  Down 

GOID 

GOTERM 

Level 

RVALUE 

c. 

Value 

GeneHits 

G0:0009628 

response  to  abiotic  stimulus 

3 

0.001886 

3.96 

7 

G0:0009653 

anatomical  structure  morphogenesis 

3 

0.00525 

2.19 

13 

G0:0043436 

oxoacid  metabolic  process 

5 

0.005573 

2.94 

8 

G0:0019752 

carboxylic  acid  metabolic  process 

6 

0.005573 

2.94 

8 

G0:0006082 

organic  acid  metabolic  process 

4 

0.006347 

2.87 

8 

G0:0042180 

cellular  ketone  metabolic  process 

4 

0.006541 

2.86 

8 

G0:0010817 

regulation  of  hormone  levels 

4 

0.007475 

4.09 

5 

G0:0048666 

neuron  development 

5 

0.008904 

2.97 

7 

G0:0048812 

neuron  projection  morphogenesis 

7 

0.009043 

3.34 

6 

G0:0006508 

proteolysis 

5 

0.014732 

2.48 

8 
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